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Abstract 

We  examine  how  to  construct  a  spatial  manifold  and  its  geometry  from  the  entanglement  struc¬ 
ture  of  an  abstract  quantum  state  in  Hilbert  space.  Given  a  decomposition  of  Hilbert  space  %  into 
a  tensor  product  of  factors,  we  consider  a  class  of  “redundancy-constrained  states”  in  %  that  gen¬ 
eralize  the  area-law  behavior  for  entanglement  entropy  usually  found  in  condensed-matter  systems 
with  gapped  local  Hamiltonians.  Using  mutual  information  to  define  a  distance  measure  on  the 
graph,  we  employ  classical  multidimensional  scaling  to  extract  the  best-fit  spatial  dimensionality  of 
the  emergent  geometry.  We  then  show  that  entanglement  perturbations  on  such  emergent  geome¬ 
tries  naturally  give  rise  to  local  modifications  of  spatial  curvature  which  obey  a  (spatial)  analog  of 
Einstein’s  equation.  The  Hilbert  space  corresponding  to  a  region  of  flat  space  is  finite-dimensional 
and  scales  as  the  volume,  though  the  entropy  (and  the  maximum  change  thereof)  scales  like  the 
area  of  the  boundary.  A  version  of  the  ER=EPR  conjecture  is  recovered,  in  that  perturbations  that 
entangle  distant  parts  of  the  emergent  geometry  generate  a  configuration  that  may  be  considered 
as  a  highly  quantum  wormhole. 


*  cjcao@caltech.edu 

t  seancarroll@gmail.com 

*  spiros@caltech.edu 


1 


CONTENTS 


I.  Introduction  [3] 

II.  Area-Law  Entanglement  0 

A.  Gravity  and  Entropy  Bounds  [5] 

B.  Area  Laws  and  Graphs  [7] 

III.  Emergent  Space  [9] 

A.  Redundancy-Constrained  States  02 

B.  Metric  from  Information  ED 

C.  Classical  Multidimensional  Scaling  [12] 

D.  Examples  with  Area-Law  States  [L4] 

IV.  Curvature  and  Entanglement  Perturbations  [T7] 

A.  Entanglement  Perturbations  [TT] 

B.  Geometric  Implications  [T8] 

1.  Effects  of  Local  Entanglement  Perturbations  [T9] 

2.  Effects  of  Nonlocal  Entanglement  Perturbations  [21] 

V.  Energy  and  Einstein’s  equation  [22 ] 

A.  Renormalization  and  the  Low  Energy  Effective  Theory  1221 

B.  Energy  and  Gravity  [24] 

VI.  Discussion  EH 

Acknowledgements  [28] 

A.  Redundancy-Constraint  and  Coarse-graining  [29 ] 

B.  Entanglement  Perturbations  and  Coarse  Curvature  1301 

References  [32] 


2 


I.  INTRODUCTION 


Quantum-mechanical  theories  are  generally  thought  of  as  theories  of  something.  Quan¬ 
tum  states  are  square-integrable  complex-valued  functions  of  the  configuration  of  some  par¬ 
ticular  kind  of  “stuff,”  where  that  stuff  may  be  a  simple  harmonic  oscillator,  a  set  of  inter¬ 
acting  spins,  or  a  collection  of  relativistic  fields. 

But  quantum  states  live  in  Hilbert  space,  a  complete  complex  vector  space  of  specified 
dimension  with  an  inner  product.  The  same  quantum  states,  even  with  the  same  dynamics, 
might  be  thought  of  as  describing  very  different  kinds  of  stuff.  Coleman  long  ago  showed 
that  the  quantum  theory  of  the  sine-Gordon  boson  in  1+1  dimensions  was  equivalent  to 
that  of  a  massive  Thirring  fermion  [Tj.  AdS/CFT  posits  an  equivalence  (in  a  certain  limit) 
between  a  conformal  field  theory  in  a  fixed  d- dimensional  Minkowski  background  and  a 
gravitational  theory  in  a  dynamical  (d+  l)-dimensional  spacetime  with  asymptotically  anti- 
de  Sitter  boundary  conditions  [2j.  The  wave  functions  of  a  single  quantum  theory  can  be 
represented  in  very  different-looking  ways.  It  is  therefore  interesting  to  consider  the  inverse 
problem  to  “quantizing”  a  theory:  starting  with  a  quantum  theory  defined  in  Hilbert  space, 
and  asking  what  it  is  a  theory  of.  In  this  paper  we  take  steps  toward  deriving  the  existence 
and  properties  of  space  itself  from  an  intrinsically  quantum  description  using  entanglement. 

A  good  deal  of  recent  work  has  addressed  the  relationship  between  quantum  entangle¬ 
ment  and  spacetime  geometry.  Much  of  the  attention  has  focused  on  holographic  models, 
especially  in  an  AdS/CFT  context.  Entanglement  in  the  boundary  theory  has  been  directly 
related  to  bulk  geometry,  including  deriving  the  bulk  Einstein  equation  from  the  entangle¬ 
ment  first  law  (EFL)  [31/6] .  (The  EFL  relates  a  perturbative  change  in  the  entropy  of  a 
density  matrix  to  the  change  in  the  expectation  value  of  its  modular  Hamiltonian,  as  dis¬ 
cussed  below.)  Tensor  networks  have  provided  a  connection  between  emergent  geometry, 
quantum  information,  and  many-body  systems  urn- 

It  is  also  possible  to  investigate  the  entanglement /geometry  connection  directly  in  a 
spacetime  bulk.  The  ER=EPR  conjecture  relates  entanglement  between  individual  particles 
to  spacetime  wormholes  P3H2U.  Consider  two  entangled  particles,  separated  by  a  long 
distance,  compared  to  the  same  particles  but  unentangled.  If  sufficient  entanglement  gives 
rise  to  a  wormhole  geometry,  some  weak  gravitational  effects  should  arise  from  small  amounts 
of  entanglement,  and  evidence  for  this  phenomenon  can  be  found  in  the  context  of  AdS/CFT 
mmm-  From  a  different  perspective,  Jacobson  has  argued  that  Einstein’s  equation  can 
be  derived  from  bulk  entanglement  under  an  assumption  of  local  thermodynamic  equilibrium 
between  infrared  and  ultraviolet  degrees  of  freedom  [23H2S]- 

While  the  current  paper  is  inspired  by  the  idea  of  emerging  space  from  entanglement,  our 
approach  of  bulk  emergent  gravity  differs  from  the  aforementioned  papers  in  that  our  starting 
point  is  directly  in  Hilbert  space,  rather  than  perturbations  around  a  boundary  theory  or 
a  semiclassical  spacetime.  We  will  first  try  to  construct  a  generic  framework  by  which  an 
approximate  sense  of  geometry  can  be  defined  purely  from  the  entanglement  structure  of 
some  special  states.  We  conjecture  that  mutual  information  (See  J2H1  27j  for  a  review), 
similar  to  suggestions  by  pm  HZ],  can  be  used  to  associate  spatial  manifolds  with  certain 
kinds  of  quantum  states.  More  tentatively,  we  explore  the  possibility  that  perturbations  of 
the  state  lead  to  relations  between  the  modular  Hamiltonian  and  the  emergent  geometry  that 
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can  be  interpreted  as  Einstein’s  equation,  as  has  been  suggested  in  a  holographic  AdS / CFT 
context.  In  doing  so  we  will  follow  some  of  the  logic  in  [5]  and  [23] -  In  particular,  we  show 
that  “nonlocal”  perturbations  that  entangle  distant  parts  of  the  emergent  geometry,  similar 
to  the  case  in  ER=EPR,  will  give  rise  to  what  might  be  understood  as  a  highly  quantum 
wormhole,  where  spatial  curvature  generated  by  (modular)  energy,  in  a  manner  similar  to 
Einstein’s  equation,  is  localized  at  the  wormhole  “mouths.” 

Our  basic  strategy  is  as  follows: 

•  Decompose  Hilbert  space  into  a  large  number  of  factors,  'H  =  'Hp.  Each  factor  is 
finite- dimensional. 

•  Consider  states  |^>0)  £  R  that  are  “redundancy-constrained,”  a  generalization  of  states 
in  which  the  entropy  of  a  region  obeys  an  area  law. 

•  Use  the  mutual  information  between  factors  A  and  B,  I(A:  B)  =  S(A)+S(B)—S(AB), 
to  define  a  metric  on  the  graph  connecting  the  factors  "Hp. 

•  Show  how  to  reconstruct  smooth,  flat  geometries  from  such  a  graph  metric  (when  it 
exists) . 

•  Consider  perturbations  |^o)  — >  IV’o)  +  and  show  these  produce  local  curvature 
proportional  to  the  local  change  in  entropy. 

•  Relate  the  change  in  entropy  to  that  in  an  effective  IR  held  theory,  and  show  how  the 
entanglement  first  law  5S  =  8(K)  (where  K  is  the  modular  Hamiltonian)  implies  a 
geometry/energy  relation  reminiscent  of  Einstein’s  equation. 


We  do  not  assume  any  particular  Hamiltonian  for  the  quantum  dynamics  of  our  state,  nor 
do  we  explore  the  emergence  of  Lorentz  invariance  or  other  features  necessary  to  claim  we 
truly  have  an  effective  quantum  theory  of  gravity,  leaving  that  for  future  work. 

We  begin  the  paper  by  reviewing  entropy  bounds  and  properties  of  entanglement  for  area- 
law  systems  in  section  [TT|  In  section  HI  we  introduce  the  notion  of  redundancy  constraint 
for  entanglement  structure  and  show  how  an  approximate  sense  of  geometry  for  can  emerge 
in  such  states.  In  particular,  we  give  a  generic  outline  of  the  procedure,  followed  by  an 
example  using  an  area-law  state  from  a  gapped  system  where  it  is  possible  to  approximately 
reconstruct  space  with  Euclidean  geometry.  Then  in  section  IV  we  discuss  the  effects  of 
entanglement  perturbations  in  terms  of  the  approximate  emergent  geometry,  and  in  |V|  show 
that  an  analog  of  the  Einstein’s  equations  can  be  derived.  Finally,  in  VI  we  conclude  with 
a  few  remarks. 

As  this  work  was  being  completed  we  became  aware  of  a  paper  with  related  goals  [28]. 
There  are  also  potential  connections  with  a  number  of  approaches  to  quantum  gravity, 
including  loop  quantum  gravity  [29],  quantum  graphity  [30],  holographic  space-time  [3T,  32], 
and  random  dynamics  [83];  we  do  not  investigate  these  directly  here. 

Throughout  this  paper,  we  will  use  d  to  denote  spacetime  dimension  and  D  =  d  —  1  for 
spatial  dimensionality. 
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II.  AREA-LAW  ENTANGLEMENT 


A.  Gravity  and  Entropy  Bounds 

The  Bekenstein-Hawking  entropy  of  a  black  hole  in  3+1  dimensions  is  proportional  to 
the  area  A  of  its  event  horizon, 

5bH  =  AG  =  27F^’  ^ 

where  we  use  h  =  c  =  1  and  the  reduced  Planck  length  is  £p  =  y/8nG.  At  a  quick  glance  this 
might  seem  like  a  surprising  result,  as  the  entropy  of  a  classical  thermodynamic  system  is  an 
extensive  quantity  that  scales  with  volume  rather  than  area.  What  does  this  imply  about 
the  Hilbert  space  describing  the  quantum  system  that  is  a  black  hole,  or  spatial  regions  more 
generally? 

Consider  a  fixed  lattice  of  qubits,  with  a  spacing  /0  and  a  linear  size  r.  The  total  number  of 
qubits  is  n  ~  ( r/£0)D ,  where  D  is  the  dimensionality  of  space,  and  the  associated  dimension 
of  Hilbert  space  is  N  =  2n.  If  the  system  is  in  a  (potentially  mixed)  state  with  density 
matrix  p,  the  von  Neumann  entropy  is  S  —  —  Tr  p  log  p.  The  maximum  entropy  of  such  a 
system  is  then  Sm ax  =  log2(iV)  =  n,  proportional  to  the  system  volume.  We  might  guess  that 
gravity  provides  an  ultraviolet  cutoff  that  acts  similarly  to  a  lattice  with  Zo  =  £p-  However, 
Bekenstein  argued  that  the  vast  majority  of  the  states  included  in  such  a  calculation  are 
physically  unattainable,  and  that  the  entropy  of  a  system  with  mass  M  and  linear  size  R  is 
bounded  by  S'  <  2ir RM  [3H  '35] .  Since  a  system  with  GM  >  R/2  undergoes  gravitational 
collapse  to  a  black  hole,  this  suggests  that  ([!])  represents  an  upper  bound  on  the  entropy  of 
any  system  in  spacetime,  a  constraint  known  as  the  holographic  bound.  If  we  were  able  to 
construct  a  higher-entropy  state  with  less  energy  than  a  black  hole,  we  could  add  energy  to 
it  and  make  it  collapse  into  a  black  hole;  but  that  would  represent  a  decrease  in  entropy, 
apparently  violating  the  Second  Law.  The  Bousso  bound  [36j  provides  a  covariant  version 
of  the  holographic  bound,  ’t  Hooft  and  Susskind  built  on  this  argument  to  suggest  the 
holographic  principle:  in  theories  with  gravity,  the  total  number  of  true  degrees  of  freedom 
inside  any  region  is  proportional  to  the  area  of  the  boundary  of  that  region  133  EE],  and 
such  a  system  can  be  described  by  a  Hilbert  space  with  dimension  of  approximately 

011111?/  ~  es  ~  \  (2) 

Meanwhile,  it  is  now  appreciated  that  area-law  behavior  for  entanglement  entropy  oc¬ 
curs  in  a  variety  of  quantum  systems,  including  many  non-gravitational  condensed- matter 
examples  [39].  Divide  space  into  a  region  A  and  its  complement  A.  A  quantum  state  \ip)  is 
said  to  obey  an  area  law  if  the  entropy  Sa  of  the  reduced  density  matrix  pa  =  Tr^  |-0)  (-0 1 
satisfies 

Sa  —  pA  +  •  •  •  ,  (3) 

where  A  is  the  area  of  the  surface  bounding  A,  and  p  is  a  constant  independent  of  A.  (Here 
and  elsewhere  in  this  paper,  entropy  equalities  should  be  interpreted  as  approximations 
valid  in  the  limit  of  large  system  size.)  This  behavior  is  generally  expected  in  low-energy 
states  of  quantum  field  theories  with  an  ultraviolet  cutoff  [30]  ill]  and  those  of  discrete 
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condensed- matter  systems  with  gapped  local  Hamiltonians  (he.,  short-range  interactions) 
[42]-  hi  conformal  field  theories,  Ryu  and  Takayanagi  showed  that  the  entanglement  entropy 
of  a  region  was  related  to  area,  not  of  the  region  itself,  but  of  an  extremal  surface  in  a  dual 
bulk  geometry  [43U45]. 

The  existence  of  an  area  law  does  not  by  itself  imply  holographic  behavior;  holography 
is  a  statement  about  the  number  of  degrees  of  freedom  in  a  region,  which  is  related  to 
the  maximum  possible  entropy,  but  not  directly  to  the  entropy  of  some  specific  state  as  in 
(J3|.  (The  AdS/CFT  correspondence  is  of  course  holographic  on  the  dual  gravity  side,  but 
the  CFT  by  itself  is  not.)  In  either  a  gapped  condensed- matter  system  or  a  QFT  with  an 
ultraviolet  cutoff  £o,  we  would  still  expect  degrees  of  freedom  to  fill  the  enclosed  volume, 
and  the  subsystem  in  A  to  have  dim  "Ha  ~  ehAo)  . 

As  the  UV  cutoff  length  is  taken  to  zero,  we  find  an  infinite- dimensional  Hilbert  space  in 
any  QFT,  and  the  entropy  of  a  region  of  space  will  generically  diverge.  Nevertheless,  QFT 
reasoning  can  be  used  to  derive  a  quantum  version  of  the  Bousso  bound  [MHIE],  by  positing 
that  the  relevant  entropy  is  not  the  full  entanglement  entropy,  but  the  vacuum-subtracted 
or  “Casini”  entropy  [49] .  Given  the  reduced  density  matrix  in  some  region  A,  and  the 
reduced  density  matrix  a  a  that  we  would  obtain  had  the  system  been  in  its  vacuum  state, 
the  Casini  entropy  is  given  by 

AS  =  S(pA )  -  S(aA)  =  —  Tr  pa  log  pa  +  Trcr^  log  aA-  (4) 

This  can  be  finite  even  when  Hilbert  space  is  infinite-dimensional  and  the  individual  entropies 
S(pa)  and  S’(cta)  are  infinite.  This  procedure  sidesteps  the  question  of  whether  the  true 
physical  Hilbert  space  is  infinite-dimensional  (and  the  holographic  entropy  bounds  refer  to 
entanglement  entropy  over  and  above  that  of  the  vacuum)  or  finite-dimensional  (and  the 
Casini  regularization  is  just  a  convenient  mathematical  trick). 

One  might  imagine  being  bold  and  conjecturing  not  only  that  there  are  a  finite  number  of 
degrees  of  freedom  in  any  finite  region,  as  holography  implies,  but  also  that  the  holographic 
bound  is  not  merely  an  upper  limit,  but  an  actual  equality  [5TTM52] .  That  is,  for  any  region  of 
spacetime,  its  associated  entanglement  entropy  obeys  an  area  law  (J3]) .  Evidence  for  this  kind 
of  area  law,  and  its  relationship  to  gravity,  comes  from  different  considerations.  Jacobson 
[53]  has  argued  that  if  UV  physics  renders  entropy  finite,  then  a  thermodynamic  argument 
implies  the  existence  of  gravity,  and  also  vice-versa.  Lloyd  [54]  has  suggested  that  if  each 
quantum  event  is  associated  with  a  Planck-scale  area  removed  from  two-dimensional  surfaces 
in  the  volume  in  which  the  event  takes  place,  then  Einstein’s  equation  must  hold. 

In  this  paper,  we  examine  quantum  states  in  a  finite-dimensional  Hilbert  space  and  look 
for  emergent  spatial  geometries,  under  the  assumptions  that  distances  are  determined  by 
mutual  information  and  that  “redundancy  constraint,”  which  reduces  to  the  usual  area-law 
relationship  of  the  basic  form  (|3]) ,  holds  when  there  exists  an  emergent  geometric  interpre¬ 
tation  of  the  state.  The  conjecture  that  arbitrary  regions  of  space  are  described  in  quantum 
gravity  by  finite-dimensional  Hilbert  spaces  represents  a  significant  departure  from  our  in¬ 
tuition  derived  from  quantum  field  theory. 

We  suggest  that  the  emergence  of  geometry  from  the  entanglement  structure  of  the  state 
can  reconcile  dim  "Ha  ~  (degrees  of  freedom  proportional  to  enclosed  volume)  with 

the  holographic  principle  in  a  simple  way:  if  we  were  to  “excite”  states  in  the  interior 
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by  entangling  them  with  exterior  degrees  of  freedom,  the  emergent  geometry  would  be 
dramatically  altered  so  that  the  system  would  no  longer  resemble  a  smooth  background 
manifold.  In  other  words,  those  degrees  of  freedom  are  only  “in  the  interior”  in  a  geometric 
sense  when  they  are  entangled  with  their  neighbors  but  not  with  distant  regions,  in  a  way 
reminiscent  of  ER=EPR. 


B.  Area  Laws  and  Graphs 

Simply  being  given  a  state  |-0O)  in  a  Hilbert  space  B  is  almost  being  given  no  information 
at  all.  Hilbert  space  has  very  little  structure,  and  we  can  always  find  a  basis  (|0n)}  for 
%  such  that  (^o|0i)  =  1  and  (t/,o|0n>i)  =  0.  To  make  progress  we  need  some  additional 
data,  such  as  the  Hamiltonian  or  a  decomposition  of  T~L  into  a  tensor  product  of  factors.  In 
this  paper  we  don’t  assume  any  particular  Hamiltonian,  but  begin  by  looking  at  states  and 
decompositions  that  give  us  a  generalization  of  area-law  behavior  for  entropy. 

To  get  our  bearings,  we  start  by  considering  systems  for  which  we  have  an  assumed 
notion  of  space  and  locality,  and  states  that  obey  an  area  law  of  the  form  (]3]) ,  and  ask  how 
such  behavior  can  be  recovered  in  a  more  general  context.  Typically  such  a  state  j^o)  is  a 
low-lying  energy  state  of  a  gapped  local  system.  Its  entanglement  structure  above  a  certain 
scale  seems  to  capture  the  space  on  which  the  Hamiltonian  is  defined  [7] .  The  entanglement 
structure  of  such  states  is  highly  constrained. 

A  remarkable  feature  of  these  states  is  that  the  entanglement  structure  above  the  said 
scale  can  be  fully  characterized  once  all  the  mutual  information  between  certain  subsystems 
are  known.  Divide  the  system  into  a  set  of  sufficiently  large  non-overlapping  regions  Ap.  We 
can  calculate  the  entropy  5(Ap)  of  each  region,  as  well  as  the  mutual  information  I(AP  :  Aq) 
between  any  two  regions.  The  system  therefore  naturally  defines  a  weighted  graph  G  = 
(V,  E),  on  which  vertices  V  are  the  regions  Ap,  and  the  edges  E  between  them  are  weighted 
by  the  mutual  information  (which  is  manifestly  symmetric).  . 

The  mutual  information  between  regions  is  a  measure  of  how  correlated  they  are.  It 
provides  a  useful  way  of  characterizing  the  “distance”  between  such  regions  because  of  its 
relation  to  correlation  functions  between  operators.  We  expect  that  in  the  ground  state  of 
a  field  theory,  correlators  of  field  operators  will  decay  as  exponentials  (for  massive  fields) 
or  power  laws  (for  massless  ones).  The  mutual  information  may  reflect  this  behavior,  as  it 
provides  an  upper  bound  on  the  correlation  function  between  two  operators.  For  the  mutual 
information  between  regions  A  and  B,  we  have: 

I(A:B)  =  S (pAB 1 1 pA  ®  pB ) 

1,  |2 

>  ^  I  Pab  ~  Pa®  Pb  \ 

^  {Tr  [{pab  ~  Pa®  Pb){OaOb)}}2 

-  2 lamp 

{(0A0B)  -  (Oa)(Ob))2 

2\\oA\\*\\oBr 


(5) 

(6) 

(7) 

(8) 


7 


FIG.  1:  An  “information  graph”  in  which  vertices  represent  factors  in  a  decomposition  of 
Hilbert  space,  and  edges  are  weighted  by  the  mutual  information  between  the  factors.  In 
redundancy-constrained  states,  the  entropy  of  a  group  of  factors  (such  as  the  shaded  region 
B  containing  "Hi  <E>  H2  <H>  'Hz  <E>  'Hi  <E>  ^5.)  can  be  calculated  by  summing  over  the  mutual 
information  of  the  cut  edges,  as  in  (|9]).  In  the  following  section  we  put  a  metric  on  graphs 


of  this  form  by  relating  the  distance  between  vertices  to  the  mutual  information,  in  (13) 


and  (14). 


We  therefore  choose  to  concentrate  on  mutual  information  as  a  way  of  characterizing  emer¬ 
gent  distance  without  picking  out  any  preferred  set  of  operators. 

Consider  grouping  a  set  of  non-overlapping  subregions  Ap  into  a  larger  region  B,  dividing 
space  into  B  and  its  complement  B.  Taking  advantage  of  the  short-ranged  entanglement 
in  such  states,  the  approximate  entanglement  entropy  of  B  can  be  calculated  using  the  cut 
function, 

S(B)  =  t  ^  I(AP:A,).  (9) 

peB.geB 

To  find  the  approximate  entanglement  of  region  B,  one  simply  cuts  all  edges  connecting  B 
and  its  complement  B.  The  entanglement  entropy  is  the  sum  over  all  the  weights  assigned 
to  the  cut  edges.  This  is  similar  to  counting  the  entanglement  entropy  by  the  bond  cutting 
in  tensor  networks,  except  in  this  special  case  where  we  are  content  with  approximate  entan¬ 
glement  entropy  for  large  regions,  a  simple  graph  representation  is  sufficient.  Comparatively, 
a  tensor  network  that  characterizes  this  state  contains  far  more  entanglement  information 
than  the  simple  connectivity  captured  by  the  graphs  considered  here. 

Our  conjecture  is  that  this  graph  information  is  enough  to  capture  the  coarse  geometry  of 
this  area-law  state.  If  we  restrict  ourselves  to  work  at  scales  for  which  S  oc  A,  all  information 
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encoded  in  the  form  of  larger-scale  entanglement  is  highly  redundant.  In  a  generic  state,  the 
mutual  information  between  all  disjoint  regions  Ap,  Aq  would  not  be  enough  to  characterize 
entanglement  entropy  of  S(ArAsAtAu )  for  r,  s,  t,  u  G  V.  Naively,  to  specify  the  entanglement 
entropy  of  all  larger  regions  B,  one  needs  on  the  order  of  0(2N)  data  points,  where  N  is  the 
number  of  vertices  (Hilbert-space  factors).  However,  in  the  special  case  of  area-law  entan¬ 
glement,  it  suffices  to  specify  all  the  mutual  information  between  N  factors.  The  amount 
of  classical  bits  needed  to  store  this  is  only  of  order  0(N2).  Therefore,  all  larger  partition 
entanglement  entropy  data  are  “redundant”  as  they  are  captured  by  the  mutual  informa¬ 
tion  of  smaller  parts.  Because  all  subsequent  higher-partition  entanglement  information  is 
encoded  in  the  mutual  information  between  all  suitably  chosen  partitions,  the  approximate 
geometric  information  above  the  chosen  scale  of  partitioning  can  be  characterized  by  the 
graph  representation. 

One  may  worry  that  the  subleading  terms  in  the  area-law  function  can  scale  as  volume 
and  therefore  ruin  the  estimate  for  higher-partition  entropies  at  some  level  of  coarse-graining. 
This  is,  however,  an  over-estimation.  The  entropy  of  a  region  with  approximate  radius  r 
computed  by  the  cut  function  assumes  a  strict  area  law,  which  scales  as  rD~l  for  a  .D-spatial- 
dimensional  area-law  system.  This  is  off  from  the  actual  entropy  by  amount  rD~ 2  +  . . . , 
where  missing  terms  have  even  lower  power  in  r.  The  relative  error,  which  scales  as  r_1, 
vanishes  in  the  large-region  limit. 

On  the  other  hand,  if  one  keeps  all  sub-leading  terms,  the  correct  edge  weights  one  should 
assign  are  given  by  the  intersecting  area  plus  an  error  term, 

I(Ap :  Hq)  =  oiA^Ap  n  Aq)  +  /3S.  (10) 

Therefore  one  may  worry  that  in  our  subsequent  estimate  of  entropy  for  a  bigger  region, 
the  error  term  may  accumulate  as  rD .  But  since  the  system  is  dominated  by  short-range 
entanglement,  the  number  of  edge  cuts  only  scales  as  rD~l.  So  in  the  worst  case  scenario, 
the  sublcading  terms  will  contribute  a  term  that  scales  as  area.  Therefore,  the  error  in  using 
the  cut  function  as  an  estimate  for  the  entropy  of  a  region  A  in  an  area-law  system  is  upper 
bounded  by  a  term  j3A  for  some  /3,  which  one  may  absorb  by  redefining  ol  —  a  +  /3.  In  this 
discussion,  we  are  not  concerning  ourselves  with  the  specific  value  of  a,  so  the  sub-leading 
terms  only  minimally  change  the  results. 


III.  EMERGENT  SPACE 

A.  Redundancy-Constrained  States 

Having  established  the  above  properties  for  area-law  states  in  systems  for  which  space 
and  locality  are  defined,  we  now  turn  to  a  more  general  context.  For  area-law  states,  the 
entanglement  information  between  different  factors  of  Hilbert  space  is  sufficiently  redundant 
that  it  can  be  effectively  characterized  by  only  limited  knowledge  of  mutual  information  [53] . 
In  the  rest  of  this  work,  we  restrict  ourselves  to  the  study  of  quantum  states  that  are  ap¬ 
proximately  “redundancy-constrained,”  defined  by  slight  generalizations  of  the  observations 
we  made  for  area-law  states  using  purely  entanglement  information. 
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Consider  a  quantum  state  |^0)  £  %  —  <S)p  Up-  We  say  that  the  state  is  redundancy- 
constrained,  (RC)  if,  for  any  subsystem  B  constructed  as  a  tensor  product  of  some  subset 
of  the  {'Hp},  its  entanglement  entropy  is  given  by  a  cut  function  of  the  form  (J9]) ,  where  Ap 
denotes  the  subsystem  that  lives  in  the  Hilbert  subspace  Up.  Note  that  there  is  no  geometric 
meaning  associated  with  the  Hilbert  space  at  this  point. 

Due  to  the  redundancy  of  the  entanglement  entropy  information,  the  entanglement  struc¬ 
ture  for  more  coarse-grained  partitions  can  be  sufficiently  captured  by  quantum  mutual  in¬ 
formation,  and  hence  admit  a  graph  description  as  in  the  previous  section.  The  vertices  of 
the  graph  label  subregions,  and  edge  weights  are  given  by  their  mutual  information.  By 
0,  it  immediately  follows  that  the  degree  of  each  vertex  Ap  (the  number  of  edges  emerging 
from  it)  is  bounded  from  above  by 

deg(Hp)  =  ^2l(Ap:Aq)  <  2 Smax(Ap)  <  2  In  Dp,  (11) 

q 

where  dim  1-La,,  =  Dp. 

RC  states  admit  the  same  graph  construction  as  area- law  states,  G  =  (V,E),  where 
vertices  are  Hilbert-space  factors  and  edges  are  weighted  by  the  mutual  information  between 
them.  Such  states  can  be  seen  as  a  straightforward  generalization  of  states  with  area-law 
scaling  that  also  lie  in  the  area-law  entropy  cone  [55j .  As  such,  they  form  a  superset  of 
area-law  states  which  also  satisfy  the  holographic  inequalities  [5BJ.  This  doesn’t  imply, 
however,  that  such  states  have  holographic  duals.  It  is  easy  to  check  that  satisfaction  of  all 
holographic  inequalities  is  not  a  sufficient  condition  to  indicate  if  a  state  has  a  holographic 
dual.  (As  a  simple  example,  we  know  that  area-law  states  from  a  gapped  system  don’t  have 
holographic  duals,  yet  they  still  satisfy  the  holographic  inequalities.) 

The  individual  Hilbert-space  factors  /HP  are  not  necessarily  qubits  or  some  other  irre¬ 
ducible  building  blocks  of  the  space.  In  particular,  they  may  be  further  factorizable,  and  are 
required  to  be  sufficiently  large  that  redundancy-constraint  becomes  a  good  approximation, 
even  if  it  would  not  hold  at  finer  scales.  In  a  phenomenologically  relevant  model,  we  would 
expect  each  factor  to  describe  not  only  the  geometry  but  the  field  content  of  a  region  of 
space  somewhat  larger  than  the  Planck  volume,  though  we  will  not  discuss  those  details 
here.  Note  that  the  RC  property  is  preserved  under  a  coarse-graining  operation  in  which  we 
decompose  Hilbert  space  into  factors  1-ip  that  are  products  of  several  of  the  original  factors 
Hp.  We  discuss  coarse-graining  more  in  appendix  [A] 

RC  states  are  highly  non-generic;  they  represent  situations  where  entanglement  is  dom¬ 
inated  by  short-range  effects.  For  example,  a  CFT  ground  state  in  D  dimensions  with  a 
holographic  dual  is  not  RC,  although  its  entanglement  data  is  still  somewhat  redundant  in 
that  one  only  needs  the  entanglement  entropy  for  balls  of  all  radii  to  reconstruct  the  AdS 
geometry  HIE].  However,  additional  data  encoded  in  the  larger  partitions  cannot  be  char¬ 
acterized  by  mere  mutual  information  between  the  partitions  Ap  for  any  coarse-graining.  In 
this  case,  the  attempt  to  define  entanglement  entropy  as  area  or  mutual  information  doesn’t 
quite  work  in  d  spatial  dimensions  any  more  because  there  is  no  simple  additive  expression 
for  S(Bx)  from  I(Ap :  Aq).  This  extra  data  for  larger  partitions  is  essential  in  constructing 
the  emergent  dimension  with  AdS  geometry. 

At  the  same  time,  if  we  have  some  dual  bulk  fields  living  in  AdS  whose  ground  state 
is  presumably  also  short  range  entangled  [ID],  then  it  may  in  turn  be  described  by  a  RC 
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state  in  AdS  with  proper  coarse-graining.  Therefore,  if  one  has  the  complete  holographic 
dictionary,  an  experiment  of  entangling  two  copies  of  CFT  to  create  a  thermofield  double 
state  has  its  dual  experiment  with  certain  constraints  in  the  bulk,  where  the  entanglement 
now  is  directly  created  in  the  bulk  and  two  copies  of  AdS  are  turned  into  a  wormhole.  Our 
general  program,  however,  does  not  rely  on  the  existence  of  a  dual  CFT. 

Although  most  states  in  ffilbert  space  given  a  certain  decomposition  are  not  redundancy- 
constrained,  RC  states  seem  like  an  appropriate  starting  point  for  investigating  theories  of 
quantum  gravity,  especially  if  area-law  behavior  for  entropy  is  universal.  For  the  remainder 
of  the  program,  we  are  going  to  focus  on  simple  RC  states  that  correspond  to  flat  space  in 
D  spatial  dimensions. 


B.  Metric  from  Information 


Consider  a  state  I'l/’o)  for  which  there  exist  a  decomposition  of  the  Hilbert  space  such 
that  l^o)  G  T-L  —  (^)p  Hp  is  redundancy-constrained.  Such  a  state  naturally  defines  a  graph 
G  =  (V,E),  with  N  vertices  labelled  by  p  and  each  edge  {p,  q}  is  weighted  by  the  mutual 
information  I(Ap:Aq).  Without  loss  of  generality,  assume  G  is  connected.  In  the  case 
when  G  has  multiple  large  disconnected  components,  one  can  simply  perform  the  procedure 
separately  for  each  connected  component. 

Our  conjecture  is  that  this  graph  contains  sufficient  information  to  define  another 
weighted  graph,  G(V,E),  on  which  the  edge  weights  can  be  interpreted  as  distances, 
thereby  defining  a  metric  space.  In  general,  passing  from  the  “information  graph”  G  to  the 
“distance  graph”  G  might  be  a  nontrivial  transformation, 

G(V,E)^G(V,E),  (12) 


with  a  different  set  of  vertices  and  edges  as  well  as  weights.  However,  we  will  make  the 
simplifying  assumptions  that  the  vertices  and  edges  remain  fixed,  so  that  the  graph  is  merely 
re- weighted,  and  furthermore  that  the  distance  weight  for  any  edge  w(p,  q )  is  determined 
solely  by  the  corresponding  mutual  information,  I(Ap  :  Aq)  (where  it  is  nonzero),  rather  than 
depending  on  the  rest  of  the  graph. 

Our  expectation  is  that  nearby  parts  of  space  have  higher  mutual  information,  while 
faraway  ones  have  lower.  We  therefore  take  as  our  ansatz  that  the  distance  between  entangled 
factors  is  some  function  of  the  mutual  information,  and  express  this  as  a  new  weight  w(p,  q) 
on  the  edges  of  our  graph.  That  is,  for  any  p,q  G  V  where  I(AP :  Aq)  ^  0,  define  the  edge 
weights  to  be 

,  ,  kc$(/(A>:A,)//o)  ip*q) 

w{p,q)=<  v  7  (13) 


(p  =  q) 


for  some  length  scale  £rc,  the  “redundancy-constraint  scale.”  No  edges  are  drawn  if 
I(Ap  :  Aq)  =  0.  Here  we  define  the  “normalized”  mutual  information  i(p  :  q)  =  I(Ap  :  Aq)/I0, 
where  normalization  Iq  is  chosen  such  that  I(Ap:  Aq) / To  =  1  when  two  regions  Ap,Aq  are 
maximally  entangled.  In  the  case  when  the  Hilbert  space  dimension  is  constant  for  all 
subregions,  we  have  I0  =  2 S(Ap)max  =  21og(dim  T-Ld)- 
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The  specific  form  of  the  scaling  function  will  presumably  be  determined  by  the  kind  of 
system  we  are  describing  ( e.g .  by  the  matter  content);  only  some  of  its  basic  properties  will 
be  crucial  to  our  considerations.  To  be  consistent  with  our  intuition,  we  require  $(1)  =  0 
and  lim.,,^0  =  oo,  namely,  the  distance  is  zero  when  two  states  are  maximally  entangled 

and  far  apart  when  they  are  unentangled.  Similar  notions  were  found  in  [10;,  [TTJ .  In  addition, 
we  choose  <3>(x)  to  be  a  non- negative  monotonically  decreasing  function  in  the  interval  [0, 1], 
where  a  smaller  mutual  information  indicates  a  larger  distance.  For  definiteness  it  may  be 
helpful  to  imagine  that  <F(x)  =  —  log(x),  as  might  be  expected  in  the  ground  state  of  a 
gapped  system  H2U57I. 

We  can  now  construct  a  metric  space  in  the  usual  way,  treating  weights  w(p,  q )  as  distances 
d(p,q).  For  vertices  connected  by  more  than  one  edge,  the  metric  d(p,q)  is  given  by  the 
shortest  distance  connecting  p  and  q.  Let  P  be  a  connected  path  between  p  and  q,  denoted 
by  the  sequence  of  vertices  P  =  (p  =  Po,Pi,P2,  ■  ■  -Pk  —  <?)■  The  metric  d(p,  q)  is  then 


k- 1 

d(p,  q)  =  min  {£  W(pn,pn+ 1)}  (14) 

n= 0 

for  all  connected  paths  P.  It  is  clear  from  the  definition  that  for  a  connected  component, 
d(p,  q)  =  d(q,p),  d(p,  q)  =  0  yy  p  =  q,  and  the  triangle  inequality  is  satisfied. 

Given  a  graph  with  N  vertices  with  a  metric  defined  on  it,  we  would  like  to  ask  whether 
it  approximates  a  smooth  manifold  of  dimension  D  <C  N.  Clearly  that  will  be  true  for  some 

graphs,  but  not  all.  One  approach  is  to  consider  an  r-ball  centered  at  p  using  the  metric  d, 

and  compute  the  entropy  of  the  reduced  density  matrix  obtained  by  tracing  out  all  regions 
outside  the  ball.  The  fractal  dimension  near  some  vertex  p  can  be  recovered  if 

S(r,p)~rDf.  (15) 


In  general  this  expression  may  not  converge  to  an  integer  Df.  In  the  case  of  integer  dimen¬ 
sion,  one  can  then  proceed  to  find  a  D  =  Df  +  1  dimensional  manifold  on  which  G  can  be 
embedded  that  comes  closest  to  preserving  the  metric  d(p,q).  We  also  assign  the  interface 
area  between  two  subregions  Ap,  Aq  as 

A:=I(Ap:Aq)/2a,  (16) 


for  some  constant  a.  Note  that  this  implies  the  area  that  encloses  the  region  Ap  in  a 
redundancy-constrained  state  is  given  by  A(Ap)  =  S(Ap)/a.  We  define  this  to  be  the 
emergent  spatial  geometry  of  the  state  and  assign  geometric  labels  to  the  Hilbert  space 
factors  based  on  the  embedding.  For  simple  geometries,  we  will  show  in  section  III C  that 
one  can  use  the  so-called  dimensionality  reduction  techniques  in  manifold  learning. 


C.  Classical  Multidimensional  Scaling 

We  now  turn  to  the  problem  of  going  from  a  graph  with  a  metric  to  a  smooth  manifold. 
One  approach  is  to  use  Regge  calculus,  which  we  investigate  in  appendix  [B}  Here  we  look  at 
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an  alternative  procedure,  multidimensional  scaling  (MDS).  For  a  more  detailed  review,  see 
e.g.  |SHJ. 


This  procedure  defines  an  embedding  of  the  graph  into  a  symmetric  manifold;  for  simplic¬ 
ity,  we  restrict  our  attention  to  cases  where  the  manifold  is  Euclidean.  The  embedding  is  an 
isometry  when  the  graph  itself  is  exactly  flat,  but  also  works  to  find  approximate  embeddings 
for  spaces  with  some  small  distortion.  In  our  current  program,  one  expects  that  there  exists 
some  natural  number  D  -C  N  where  the  corresponding  embedding  in  D-Euclidean  space  is 
(approximately)  isometric,  but  there  can  be  distortion  since  there  is  some  arbitrariness  in 
our  choice  of  the  distance  function  <f>  appearing  in  (13). 

Consider  the  distance  graph  G  =  (V,E),  with  edges  weighted  by  the  metric  distance 
d(p,q).  These  vertices  and  distances  now  define  a  metric  space  (V,d).  The  first  thing  we 
can  do  is  define  the  emergent  dimension  of  this  discrete  space.  Consider  a  subset  of  vertices, 
X  =  {n0,Ui, . . . ,vr}  C  G,  equipped  with  its  induced  metric.  X  is  a  metric  subspace  and  a 
r-simplex  of  V.  Now  construct  the  matrix 


Rij  =  ~(d(vi,v0)2  +  d(vj,v o)2  -  d(vi,Vj)2). 


(17) 


Since  the  determinant  det(R)  =  R(v 0,th, . . .  ,vr)  is  a  symmetric  function,  define  simplicial 
volume 

volJX )  =  —Wdet(R),  (18) 

T ! 

which  is  nothing  but  the  spatial  volume  of  the  r-simplex  if  X  is  a  subset  of  Euclidean  space 
equipped  with  the  induced  Euclidean  metric.  The  dimension  of  a  metric  space,  if  it  exists, 
is  the  largest  natural  number  k  for  which  there  exists  a  D-simplex  with  positive  volume.  As 
demonstrated  by  [59],  the  metric  space  can  be  isometrically  embedded  into  Euclidean  space 
with  dimension  d  if  and  only  if  the  metric  space  is  flat  and  has  dimension  <  D. 

The  output  of  MDS  applied  to  N  vertices  with  distances  d(p,  q )  embedded  into  a  D- 
dimensional  space  is  an  IV  x  .D  matrix  X,  which  can  be  thought  of  as  the  embedded  coor¬ 
dinates  of  all  the  vertices:  the  nth  row  contains  the  D  coordinate  values  of  the  nth  vertex, 
up  to  isometric  transformations. 

To  see  how  this  might  work,  imagine  for  the  moment  working  backwards:  given  some 
coordinate  matrix  X,  how  is  it  related  to  the  distances  d(p,  q)l  First  define  an  n  x  n  matrix 
B  =  XX*  =  (XO)(XO)*,  which  is  equivalent  for  coordinate  matrix  X  up  to  some  arbitrary 
orthonormal  transformation  O.  Then  we  notice  that  the  Euclidean  distances  between  two 
rows  of  X  can  be  written  as 


d{p,  q)2  =  -  xqry 

r= 1 
d 

^  ^  [XpyXpr  “f"  XqrXqr  2 XprX( 


r= 1 


Bpp  Bqq  2Bpq. 


(19) 

(20) 
(21) 


Therefore,  if  B  can  be  recovered  only  from  the  Euclidean  distances  d(p,  q ),  a  solution  for  X 
can  be  obtained. 
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The  solution  X  for  the  embedding  coordinates  is  non-unique  up  to  isometric  transfor¬ 
mations.  To  get  a  unique  solution,  we  first  impose  the  following  constraints  such  that  the 
embedding  is  centered  at  the  origin, 

N 

J2xPr  =  0,  Vr.  (22) 

p= i 

Then  it  follows  that  Bpq  =  0  and 

N  N  N 

Bpi  =  ~\  fa,  q)2  -  ~  d(p,  l)2  -  m )2)  •  (23) 

1=1  1=1  l,m=  1 

This  defines  the  components  Bpq  in  terms  of  the  graph  distances. 

We  diagonalize  B  via 

B  =  VAV*.  (24) 

Here  A  is  a  diagonal  matrix  with  diagonal  eigenvalues  Ai  >  A2  >  . . .  A n  arranged  in  descend¬ 
ing  order.  In  addition,  because  B  has  rank  D,  we  choose  the  D  eigenvectors  corresponding 
to  the  D  non-zero  eigenvalues.  A  solution  for  X  is  then 

X=(V^vi,...,vW  (25) 


which  is  an  isometric  embedding  of  N  points  into  a  D- dimensional  Euclidean  space. 

For  the  case  where  exact  embedding  is  not  possible,  i.e.,  the  distance  function  is  Eu¬ 
clidean  but  with  some  small  deviations,  there  will  be  D  dominant  positive  eigenvalues  fol¬ 
lowed  by  smaller  non-zero  eigenvalues.  We  consider  the  D-dimensional  embedding  to  be 
approximately  valid  if  =  1  —  Y2f=i  |Ai|/Eili  |Aj|  is  sufficiently  small.  One  can  quantify 
the  distortion  from  exact  embedding  in  various  ways.  For  instance,  for  the  classical  MDS 
algorithm  we  use  here,  a  so-called  stress  function  is  used  as  a  measure  of  distortion 


Stress 


(EP,g  d(p,g)dE(xp,xq ))2 
d(p,  q)2  J2P,q  dE(xp,  xq )2  ’ 


(26) 


where  dE(xp,  xq)  is  the  Euclidean  distance  between  the  two  corresponding  vertices  p,  q  in  the 
embedded  space.  Essentially,  MDS  analytically  generates  a  set  of  embedding  coordinates  in 
a  lower  dimensional  Euclidean  space  [58 >  IBP] ,  where  the  algorithm  seeks  an  optimal  Euclidean 
embedding  such  that  the  inter-vertex  distances  are  best  preserved  in  sense  that  Stress  is 
minimized.  Although  it  also  works  well  for  graph  embedding  in  highly  symmetric  surfaces 
(hyperbolic  and  spherical  as  well  as  flat)  [60l  [62] .  it  is  considered  a  difficult  problem  to  find 
an  embedding  for  generic  curved  manifolds.  The  matching  problem,  although  non-trivial, 
can  be  significantly  simplified  if  the  embedding  manifold  is  known  [63]. 


D.  Examples  with  Area-Law  States 

Let’s  see  how  MDS  works  in  practice  for  our  redundancy-constrained  quantum  states.  To 
do  this,  we  will  examine  states  whose  geometric  interpretation  is  known,  and  show  that  our 
procedure  can  recover  that  geometry. 
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FIG.  2:  Multidimensional  scaling  results  for  the  1  -d  antiferromagnetic  Heisenberg  chain. 


On  the  left  we  plot  the  eigenvalues  of  the  matrix  (24).  The  fact  that  the  first  eigenvalue  is 
much  greater  than  the  others  indicates  that  we  have  a  1  -d  embedding.  On  the  right  we 
show  the  reconstructed  geometry  by  plotting  the  first  three  coordinates  of  the  graph 

vertices. 


We  start  by  imagining  that  an  unsuspecting  group  of  theorists  have  been  handed  the 
state  l^o)  6  H  —  0^  Tip,  which  is  actually  the  ground  state  of  a  d-dimensional  gapped 
local  Hamiltonian  that  also  satisfies  an  area  law.  Although  the  theorists  are  only  given  \i/jq) 
and  its  Hilbert  space  decomposition,  they  are  tasked  with  finding  an  approximate  geometry 
for  the  state. 

We  start  by  constructing  the  graph  G  =  (V,  E ),  where  the  vertices  are  labelled  by  subre¬ 
gions  Ap,  and  edge  weights  are  given  by  the  distance  function  w(p,  q )  =  £rc$[I(Ap  :  Aq)/I0], 
as  in  (13).  For  convenience  we  choose  3>(x)  =  —  ln(x).  This  function  is  chosen  as  most 
finitely  correlated  states  have  fast  decaying  correlation  which,  in  the  limit  of  large  distances, 
is  exponentially  suppressed  m  In  particular,  this  is  satisfied  for  any  system  with  a  spectral 
gap  whose  observables  commute  at  large  distances  1523-  The  correlation  of  any  state  that 
is  locally  entangled  (finitely  correlated  states),  e.g.,  ones  that  can  be  expressed  in  terms  of 
MPS  or  PEPS  tensor  networks,  are  expected  to  take  this  form. 

Two  examples  are  illustrated  here,  corresponding  to  a  state  living  on  a  one-dimensional 
line  and  one  living  on  a  two-dimensional  plane.  In  both  cases  we  start  with  a  known  vacuum 
state  of  a  gapped  local  system,  where  correlations  are  expected  be  short-ranged.  Computing 
the  mutual  information  for  a  quantum  state  of  such  systems  is  in  general  not  an  easy  task. 
Consequently,  we  did  not  calculate  directly  the  mutual  information  from  density  matrices 
in  the  1  -d  case  and  instead  used  the  correlation  function  as  a  proxy. 

Our  one- dimensional  example  is  the  ground  state  of  an  (S'=l)  1-dimensional  antiferro¬ 
magnetic  Heisenberg  chain  m  165] .  Recalling  (J8]) ,  we  use  the  magnitude  squared  of  the 
correlation  function  as  an  estimate  for  the  mutual  information.  The  ground  state  correla- 
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FIG.  3:  Multidimensional  scaling  results  for  a  coarse-grained  2-d  toric  code  ground  state. 


Again,  the  left  shows  the  eigenvalues  of  (24)  and  the  right  shows  the  reconstructed 
geometry.  The  two  dominant  eigenvalues  show  that  the  geometry  is  two-dimensional, 
though  the  fit  is  not  as  close  as  it  was  in  the  1  -d  example.  Similarly,  the  reconstructed 

geometry  shows  a  bit  more  distortion. 


tion  function  |(S'“S'“)|  is  approximately  proportional  to  the  modified  Bessel  function  K0(r/£), 
where  a  =  x,y,z.  This  is  a  fast-decaying  correlation,  scaling  as  as  exp(— r/£)/y/r/£  in  the 
asymptotical  limit  when  r  >  (.  (For  a  =  z  the  correlator  is  supplemented  with  an  extra 
term  of  the  same  order,  given  by  ~  2£A'0(r/2£)A\(r/2£)/r.  The  distortion  is  still  minimal, 
and  in  fact  yields  a  slightly  better  isometric  embedding.) 

We  constructed  a  graph  of  100  vertices  and  assigned  edge  weights  given  by  the  square 
of  correlator.  No  coarse-graining  is  performed.  Applying  MDS  to  this  graph  returns  a 
vector  of  embedding  coordinates  in  Euclidean  space.  As  expected,  there  is  distortion  in  the 
embedding  and  the  coordinate  matrix  X  has  rank  greater  than  one.  However,  the  distortion, 
as  measured  using  eigenvalues,  has  e\  =  0.0167,  from  which  we  determine  that  it  only  slightly 
deviates  from  an  embedding  in  Id.  In  figure  [2j  a  patch  of  approximately  50  points  is  plotted. 
The  one-dimensional  nature  of  the  reconstructed  geometry  is  evident. 

Our  second  example  reconstructs  a  patch  of  the  2-d  toric  code  [66] ,  where  it  is  possible  to 
exactly  calculate  the  entropy  for  different  subregions  isa.  In  a  coarse-graining  where  each 
region  is  homeomorphic  to  a  plaquette,  the  exact  entropy  is  S  =  £  ab  —  1  =  Lqa  —  n2  —  2 n3  — 
1  =  ni  +  n2  +  n3  —  1,  where  LqA  is  the  length  of  the  boundary  that  separates  bipartitions  A 
and  B.  Here,  rq  denotes  the  number  of  sites/star  operators  that  have  i  nearest  neighbors  in 
A.  The  mutual  information  used  for  the  network  is  again  given  by  the  length  of  overlapping 
boundary,  up  to  a  constant  correction  term.  As  neighboring  spins  are  uncorrelated,  the 
constant  entanglement  entropy  offset  only  changes  the  overall  definition  of  length  scale  by  a 
constant  factor  for  a  near-uniform  coarse-graining.  As  a  result,  the  geometric  reconstruction 
of  a  2-d  patch  is  given  by  ([3]),  up  to  coordinate  rescaling.  The  distortion  now  is  visibly  higher 
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because  —  ln(x)  is  no  longer  an  ideal  ansatz  for  $(x).  The  distortion  factor  as  measured  by 
eigenvalues  for  an  embedding  in  2-d  has  62  =  0.41,  but  the  2-d  nature  of  the  emergent  space 
is  evident  from  figure  [3j 


IV.  CURVATURE  AND  ENTANGLEMENT  PERTURBATIONS 

In  this  section  we  examine  the  effects  on  our  reconstructed  spatial  geometries  of  per¬ 
turbing  the  entanglement  structure  of  our  states.  As  we  are  only  considering  space  rather 
than  spacetime,  we  cannot  directly  make  contact  with  general  relativity;  in  particular,  we 
can  say  nothing  about  the  emergence  of  dynamical  fields  obeying  local  Lorentz  invariance. 
Nevertheless,  we  will  see  that  the  induced  geometry  responds  to  perturbations  in  a  way 
reminiscent  of  Einstein’s  equation,  suggesting  that  an  emergent  spacetime  geometry  could 
naturally  recover  gravity  in  the  infrared. 


A.  Entanglement  Perturbations 

Consider  some  unperturbed  “vacuum”  density  operator  a  =  |'0o)(V’o|  £  L(jH),  for  which 
there  exists  a  D-dimensional  geometric  reconstruction  as  discussed  in  the  last  section.  (Here 
L{%)  denotes  the  space  of  complex-valued  linear  operators  on  'H,  of  which  the  density  opera¬ 
tor  is  an  element.)  We  choose  a  vertex  p  on  the  distance- weighted  graph  G  that  is  associated 
with  some  subregion  Ap  of  the  emergent  geometry.  The  reduced  density  matrix  associated 
with  the  region  is  defined  in  the  usual  way:  aAp  —  Trj  [c],  where  Ap  is  the  complement  of 
Ap.  The  entropy  of  such  a  region  is  again  SAp(crAp)  =  -  Tr [aAp  logaAJ  =  1/2  I(Ap  :  Ap). 
The  interface  area  between  regions  Ap  and  Aq  is  defined  as  a  A  =  I(Ap:Aq)/  2,  and  the 
distance  between  vertices  p,  q  is  defined  by  d(p,  q )  =  lp$(i(Ap  :  Aq)).  Recall  that  the  normal¬ 
ized  mutual  information  is  i(p:q)  =  / (Ap  :  Aq) / I0(p :  q) ,  where  i(p:q )  =  1  when  subsystems 
Ap,  Aq  are  maximally  entangled. 

There  are  a  variety  of  entanglement  perturbations  one  can  consider.  We  can  separately 
investigate  “local”  perturbations  that  change  the  entanglement  between  Ap  and  nearby 
degrees  of  freedom,  and  “nonlocal”  ones  that  introduce  entanglement  between  Ap  and  degrees 
of  freedom  far  away;  the  latter  can  be  modeled  by  nonunitary  transformations  on  'Ha  . 

A  local  perturbation  is  generated  by  some  unitary  operator  UAA  acting  on  the  original 
system  T~L  =  HAp  <S>  HAp.  The  perturbed  state  is  p  =  A  crUA  A  .  From  the  definition  of 
mutual  information,  we  know  that 

5I(AP :  Ap)  =  5SAp  +  5SAp  -  SSApAp 

=  2  5SAp,  (local)  (27) 

where  6 Si  =  Si(p)  —  Si(cr)  denotes  the  infinitesimal  change  of  entanglement  entropy  for  region 
i.  (This  relation  also  holds  for  finite  changes  in  entropy.)  The  second  equality  follows  because 
SSAp  =  SSAp  and  SSApAp  =  0,  since  UApAp  does  not  change  the  total  entropy  of  the  system. 
By  construction,  the  definitions  of  area  and  length  are  related  to  mutual  information  of  the 
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quantum  state;  as  we  will  soon  discover,  the  entanglement  perturbation  here  is  tantamount 
to  a  local  curvature  perturbation  at  Ap. 

Nonlocal  entanglement  perturbations  correspond  to  applying  a  lossy  quantum  channel 
A,  which  can  equivalently  be  treated  as  a  completely  positive  and  trace  preserving  (CPTP) 
map,  to  the  system.  To  that  end  we  introduce  an  extended  Hilbert  space 


U*  ='Hap®'Hap®'Ub,  (28) 

where  B  represents  some  ancillary  degrees  of  freedom  that  are  initially  unentangled  with 
those  in  Ap.  One  can  think  of  B ,  described  by  state  aB  e  L(W.b),  as  a  different  patch  of 
emergent  space,  or  simply  some  degrees  of  freedom  that  the  system  has  not  yet  encountered. 
A  nonlocal  perturbation  is  enacted  by  a  unitary  UApApB  =  Uapb  <S>  IAp  that  acts  only  on 
the  degrees  of  freedom  in  Ap  and  B.  The  perturbed  state  is  p  =  A(cr)  =  Tyb[U\  .  a  ® 

a BUApApB\-  la  this  case,  the  change  in  mutual  information  between  Ap  and  Ap  is  non¬ 
positive,  and  will  depend  on  the  local  entanglement  structure  as  well  as  the  entangling 
unitary. 

Let  F\(ASap]  Ap,  Ap)  be  a  function  that  describes  the  finite  change  in  mutual  information 
between  Ap  and  its  complement.  The  specific  implementation  of  this  function  will  depend 
on  A  and  the  entanglement  structures  related  to  the  regions  of  interest,  Ap  and  Ap.  The 
total  change  in  mutual  information  in  this  case  is  given  by 


AI(AP  :  Ap)  —  F\(ASap]  Ap,  Ap). 


(29) 


Because  the  change  in  mutual  information  has  to  be  zero  when  no  unitary  is  applied,  we 
must  have  Pa(0;  Ap,  Ap)  =  0.  For  infinitesimal  perturbations,  we  can  write 


SI(AP :  Ap) 


5S 


dF\(ASAp ;  Ap,  Ap) 


d(ASAp) 
SSApF^0;Ap,Ap). 


ASap=0 


(nonlocal) 


(30) 

(31) 


We  have  written  the  change  in  mutual  information  as  if  it  is  proportional  to  the  change  in 
entropy,  but  note  that  (in  contrast  with  the  case  of  local  perturbations)  here  the  propor¬ 
tionality  is  not  a  universal  constant,  but  rather  a  factor  that  depends  on  the  channel  A. 
In  general,  for  mixed  states  a,  aB  which  are  not  maximally  entangled,  we  can  easily  find 
unitary  operations  where  FA(0]  Ap,  Ap)  ^  0.  Note  that  this  relation  differs  from  the  local 
perturbation  by  an  F^-dependent  constant  factor. 


B.  Geometric  Implications 

We  now  consider  the  effect  of  an  entanglement  perturbation  on  the  emergent  spatial  ge¬ 
ometry.  In  this  section  we  imagine  mapping  our  graph  to  a  Riemannian  embedding  manifold, 
as  we  did  for  the  vacuum  case  using  MDS  in  the  previous  section.  In  appendix  [B]  we  study 
the  problem  using  Regge  calculus. 
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FIG.  4:  For  a  graph  G  embedded  in  some  manifold,  we  assign  subregion  Ap  (dark  blue 
region)  to  the  vertex  p,  which  is  connected  to  adjacent  vertices  q  (Black  solid  line).  An 
entanglement  perturbation  that  decreases  the  mutual  information  between  Ap  with  its 
neighbors  elongates  the  connected  edges  (dashed  red  lines),  creating  an  angular  deficit 
which  is  related  to  the  curvature  perturbation  at  p. 


Although  it  is  operationally  difficult  to  find  for  the  graph  an  embedding  manifold  with 
variable  curvature,  it  is  considerably  more  tractable  if  we  only  wish  to  quantify  a  perturba¬ 
tion  around  a  known  manifold  that  corresponds  to  the  density  matrix  a.  Namely,  in  order 
for  the  manifold  to  be  a  good  embedding,  its  perturbed  form  should  at  least  be  consistent 
with  the  deviations  in  area  and  geodesic  lengths.  Since  we  have  outlined  explicit  algorithms 
for  flat  space  configurations,  here  we  assume  that  a  D- dimensional  flat  configuration  A4  has 
been  obtained  using  the  above  framework. 


1.  Effects  of  Local  Entanglement  Perturbations 

We  begin  by  considering  a  local  perturbation  that  decreases  the  entropy  of  our  region, 

SSa,  <  0,  (32) 

which  as  we  will  see  induces  positive  spatial  curvature.  Thus  we  are  considering  a  local 
operation  that  decreases  the  entanglement  between  Ap  and  the  rest  of  Hilbert  space.  Without 
altering  the  dimension,  the  minimal  change  to  A4  that  can  be  imposed  is  some  perturbation 
to  spatial  curvature  at  p.  For  the  simplest  case,  let  Ap  be  a  region  that  contains  a  single 
graph  vertex  p  whose  entanglement  with  the  adjacent  vertices  q  (regions  Aq)  is  gradually 
decreased. 

Following  Jacobson  [23],  we  proceed  by  defining  Riemann  normal  coordinates  in  the 
vicinity  of  Ap, 

hij  =  8i:j  -  ^r2Rijklxkxl  +  0(r3).  (33) 

Consider  the  perturbed  subregion  fixed  at  some  constant  volume  V,  of  characteristic  linear 
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size  r  =  VX!D .  The  decrease  in  area  under  the  perturbation  is  given  by 


SA 


QD_irD+1 

2D(D  +  2 ) 


7? 

/VP? 


(34) 


where  1ZP  =  R,tj13  (p)  is  the  spatial  curvature  scalar  and  Q^-i  is  the  volume  of  a  unit  ( D  —  1)- 
spliere.  We  know  that  the  boundary  area  is  defined  by  local  mutual  information,  namely, 

SI(Ar:A,)  »  (35) 

Aq€.Ap 


Because  the  system  is  only  short-range  entangled,  the  local  mutual  information  is  well- 
approximated  by  the  mutual  information  between  the  region  and  its  complement.  For  in¬ 
stance,  if  the  graph  that  captures  entanglement  structure  for  the  toric  code  ground  state 
above  the  RC  scale  is  used,  then  the  two  quantities  will  be  exactly  equal.  In  general  for  sys¬ 
tems  with  exponentially  decaying  mutual  information,  the  error  with  this  estimation  is  also 
upper-bounded  by  a  quantity  of  order  exp(— t/Irc),  which  vanishes  as  long  as  the  vertices 
correspond  to  sufficiently  coarse-grained  regions. 

Plugging  (27)  for  the  change  in  mutual  information  due  to  an  infinitesimal  local  pertur¬ 
bation  into  (34)  and  (35),  we  can  relate  the  curvature  scalar  to  the  entropy  perturbation 
by 

2D(D  +  2) 


7lp  = 


\d+i 



(36) 


O' fin- i(7^rcJ 

Here  we  have  set  r  equal  to  qf'Rc,  a  characteristic  size  of  the  region  for  some  constant  7. 
Such  an  approximation  is  most  accurate  when  the  symmetry  is  also  approximately  reflected 
by  the  graph.  Because  SSap  <  0,  the  induced  curvature  is  positive. 

The  relation  between  curvature  and  the  entropy  perturbation  can  alternatively  be  derived 
by  using  (35)  to  estimate  the  decrease  in  mutual  information  for  each  individual  edge,  and 
relate  that  to  a  length  excess.  Since  the  edge  weights  sum  to  the  total  change  in  entropy, 
for  each  edge  we  can  write 

5I(Ap:Aq)  =  2Vq5SAp,  (37) 


for  some  constants  r]q  such  that  ^2qr]q  =  1.  The  values  of  r]q  are  determined  by  the  graph 
structure  near  p ,  as  well  as  the  unitary  UAA.  For  a  unitary  that  symmetrically  disentangles 
all  the  edges  on  a  regular  lattice,  r)q  —  1/  deg(p). 

Alternatively,  we  can  also  relate  the  change  in  entropy  to  the  curvature  perturbation  by 
considering  the  change  in  linear  size.  The  radius  excess  for  the  same  perturbation  at  some 
fixed  area  is 


5d 


r3Tlp 

6D(D  +  2)  ’ 


(38) 


Recall  that  distance  is  related  to  mutual  information  by  d(p,q )  =  £-Rc&(i(p  ■  <?)),  where  the 
normalized  mutual  information  is  i(p :  q)  =  I(Ap  :  Aq)/I0(p  :  q).  Assuming  an  approximately 
symmetric  configuration,  to  leading  order  we  have 


Sd(p,  q)  =  Trc& (i(p :  q))6i(p :  q)  +  0(5i2)  (39) 

=  -2riqTRc\&{i(p\q))\5SA/h{p-q)  +  0(6i2), 
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where  for  the  last  line  we  used  (27)  to  relate  the  linear  change  in  entropy  to  the  change  in 
the  distance  function.  Note  that  <f>'  =  d^/di  is  negative  by  construction. 

Comparing  (38)  to  (39)  yields  an  alternative  relation  between  curvature  and  entanglement 
entropy, 

l2r)q\$'(i(p:q))\D{D  +  2) 

- OJAp, 


Tip  = 


To(p-qyid3 


(40) 


where  again  we  have  set  r  =  y^Rc-  One  can  check  that  for  1/a  oc  >  the  fwo  results  (36) 
and  (40)  are  equivalent  up  to  some  dimension-dependent  choice  of  r)q,  the  inverse  function  $, 
and  constant  factor  a.  Both  imply  positive  curvature  for  local  disentangling  perturbations. 
Similarly,  entangling  perturbations  with  5SV  >  0  yields  negative  spatial  curvature. 


2.  Effects  of  Nonlocal  Entanglement  Perturbations 


The  derivation  with  nonlocal  entanglement  perturbation  is  similar,  where  we  simply  re¬ 
place  the  constant  proportionality  factors  with  the  channel-dependent  factor  0).  Repeat¬ 
ing  the  above  steps,  analogously  to  the  area  deficit  condition  ( 36 )  we  have 

D(D  +  2)F^(0;  Ap,  Ap) 


Tip  = 


olVLd- i(7^rc)d+1 


-5S 


Ap  i 


(41) 


while  analogously  to  the  radius  deficit  (40)  we  obtain 


Tip  = 


6Ff(0-1Ap,Aq)\V(i(p:q))\D(D  +  2) 
To(p  ■  <?Kr c73 


SSjl,. 


(42) 


Interestingly,  we  find  that  nonlocal  entanglement  perturbations  are  only  able  to  generate 
positive  curvature  perturbations.  Because  SI(Ap  :  Ap)  <  0  under  any  operations  acting  on 
Ap  and  B ,  from  (31)  we  must  have 


*a(0;  Ap,  Ap)  <  0 


(43) 


for  a  generic  entangling  unitary  when  SSa  >  0.  If  SSap  <  0,  then  it  must  follow  that 
F\( 0;  Ap,  Ap)  >  0  for  the  same  reason. 

The  nonlocal  case  is  also  interesting  due  to  its  connection  with  the  ER=EPR  conjecture. 
In  this  case,  some  spatial  region  Ap,  described  by  some  mixed  state  gav ,  is  far  separated 
from  some  other  spatial  region  B ,  with  corresponding  mixed  state  <Jr.  An  entangling  unitary 
then  creates  some  weak  entanglement  between  the  regions,  similar  to  having  an  EPR  pair 
shared  between  them.  Such  entanglement,  as  we  have  seen,  decreases  the  mutual  information 
between  Ap  (B)  and  their  respective  neighboring  regions.  This,  we  claim,  can  be  interpreted 
as  a  quantum  proto-wormhole.  No  smooth  classical  geometry  is  present  to  form  the  usual 
ER  bridge;  nevertheless,  the  entanglement  backreacts  on  the  emergent  geometry  in  a  way 
such  that  positive  modular  energy  “curves”  the  spatial  regions  near  the  “wormhole  mouths.” 

A  large  entanglement  modification  beyond  the  perturbative  limit  does  not  create  discon¬ 
nected  regions.  Heuristically,  because  the  entanglement  of  Ap  is  always  constant  under  such 
a  unitary,  it  must  become  more  entangled  with  region  B.  From  the  point  of  view  of  emergent 
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FIG.  5:  For  perturbations  that  slightly  entangle  two  regions  of  the  emergent  space,  as 
represented  by  the  vertices,  a  positive  curvature  perturbation  is  induced  locally  near  each 
perturbed  site.  We  may  interpret  this  as  a  highly  “quantum  wormhole.”  The  dotted  red 
line  joining  p  and  p'  denotes  some  trace  amount  of  entanglement  between  the  two 

subsystems. 


geometry  from  entanglement,  it  implies  that  the  region  Ap  should  also  be  connected  to  the 
distant  region  B  in  some  way.  When  such  a  connection  between  the  two  regions  becomes 
manifestly  geometric,  the  process  may  then  be  interpreted  as  the  formation  of  a  classical 
wormhole. 

However,  because  the  function  F\  depends  on  the  entanglement  structure  of  pAA  and  A 
in  addition  to  A Sa,  the  entropy-curvature  relation  does  not  seem  to  be  universal  as  in  the 
local  case.  Interpreting  nonlocal  effects  as  “gravitational”  in  this  model  may  be  in  tension 
with  our  expectations  for  a  theory  of  gravity,  although  further  assumptions  on  symmetries 
in  entanglement  structure  may  resolve  this  issue. 


V.  ENERGY  AND  EINSTEIN’S  EQUATION 

We  have  seen  how  spatial  geometry  can  emerge  from  the  entanglement  structure  of  a 
quantum  state,  and  how  that  geometry  changes  under  perturbations.  This  is  a  long  way 
from  completely  recovering  the  curved  spacetime  of  general  relativity,  both  because  we  don’t 
have  a  covariant  theory  with  dynamics,  and  because  we  haven’t  related  features  of  the  state 
to  an  effective  stress-energy.  We  can  address  the  second  of  these  points  by  considering 
general  features  of  a  map  from  our  original  theory  to  that  of  an  effective  held  theory  on  a 
fixed  spacetime  background,  then  appealing  to  the  entanglement  first  law  (EFL);  we  leave 
the  issue  of  dynamics  to  future  work.  Our  approach  here  is  similar  in  spirit  to  previous  work 
in  AdS/CFT  [3H6]  and  bulk  entropic  gravity  [23b25]. 


A.  Renormalization  and  the  Low  Energy  Effective  Theory 

To  make  contact  with  semiclassical  gravitation,  we  need  to  understand  how  an  effective 
held  theory,  and  in  particular  a  local  energy  density,  can  emerge  from  our  Hilbert-space 
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formalism.  This  is  a  nontrivial  problem,  and  here  we  simply  sketch  some  steps  toward 
a  solution,  by  integrating  out  ultraviolet  (gravitational)  degrees  of  freedom  to  obtain  an 
infrared  field  theory  propagating  on  a  background.  We  argue  that,  for  local  entanglement 
perturbations,  a  perturbation  that  decreases  the  entropy  of  a  quantum-gravity  state  in  a 
region  A  will  correspond  to  increased  entropy  density  in  the  effective  1R  held  theory,  which 
lives  in  a  lower-dimensional  (coarse-grained)  Hilbert  space. 

Our  construction  posits  an  RC  scale  £RC  at  distances  greater  than  which  the  state  obeys 
the  redundancy-constraint  condition.  Intuitively  we  expect  £rc  to  be  close  to  the  reduced 
Planck  length  tp,  at  which  the  spatial  geometry  has  barely  emerged.  Let  the  corresponding 
UV  energy  scale  be  Arc  =  1/Uic-  We  imagine  an  “emergence”  map 

S  :  H  -)•  Heft(Arc),  (44) 


which  maps  states  in  the  Hilbert  space  of  our  original  theory  to  those  of  an  effective  held 
theory  with  cutoff  Arc  in  a  semiclassical  spacetime  background. 

To  study  the  relationship  of  entropy  and  energy  in  this  emergent  low-energy  effective 
description,  we  consider  the  RG  how  of  the  theory,  defined  by  a  parameterized  map  that 
takes  the  theory  at  Arc  and  hows  it  to  a  lower  scale  A  by  integrating  out  UV  degrees  of 
freedom: 

Tx  :  Heft(Arc)  -)•  %eft(a).  (45) 


This  how  is  defined  purely  in  the  context  of  QFT  in  curved  spacetime,  so  that  the  background 
geometry  remains  hxed.  Note  that  we  could  also  discuss  RG  how  directly  in  the  entanglement 
language,  where  it  would  be  enacted  by  a  tensor  network  similar  to  MERA  [68] .  There,  the 
equivalent  of  JR  would  be  a  quantum  channel  that  could  be  defined  by  a  unitary  circuit  if 
we  include  ancillae  representing  UV  degrees  of  freedom  that  are  integrated  out.  In  this  case, 
A  depends  on  the  number  of  layers  of  the  MERA  tensor  network. 

Given  that  (36)  relates  a  spatial  curvature  perturbation  in  the  emergent  geometry  to  a 
change  in  entropy  in  the  full  theory,  we  would  like  to  know  how  this  entropy  change  relates 
to  that  of  the  vacuum-subtracted  entropy  JS'EFT(V  in  the  effective  held  theory  with  cutoff 
A  defined  on  a  background  (and  ultimately  to  the  emergent  mass/energy  in  that  theory). 
We  posit  that  they  are  related  by  a  positive  constant  K\  that  depends  on  the  cutoff  but  is 
otherwise  universal: 

SS  =  -kx5Seftw.  (46) 


The  minus  sign  deserves  some  comment.  A  perturbation  that  disentangles  a  Hilbert 
space  factor  Ap  decreases  its  entropy,  SSap  <  0,  while  inducing  positive  spatial  curvature  in 
the  emergent  geometry  and  a  decrease  in  the  area  Ap  of  the  boundary  of  the  corresponding 
region.  Naively,  a  decrease  in  boundary  area  results  in  the  decrease  of  the  entropy  of  a  region 
in  a  cutoff  effective  held  theory.  However,  that  expectation  conies  from  changing  the  area 
of  a  near-vacuum  state  in  a  hxed  background  geometry.  Here  we  have  a  different  situation, 
where  the  perturbation  affects  both  the  geometry  and  the  EFT  state  defined  on  it. 

In  that  context,  as  Jacobson  has  argued  [23],  we  expect  an  equilibrium  condition  for 
entanglement  in  small  regions  of  spacetime.  Consider  a  perturbation  of  the  EFT  defined  on 
a  semiclassical  background,  which  changes  both  the  background  geometry  and  the  quantum 
state  of  the  helds.  Fix  a  region  in  which  we  keep  the  spatial  volume  constant.  In  the 
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spirit  of  holography  and  the  Generalized  Second  Law,  the  total  entropy  in  the  region  can  be 
considered  as  the  sum  of  an  area  term  representing  UV  quantum  gravity  modes  plus  a  term 
for  the  IR  effective  field  theory, 

*Stotal  r]dvA  +  SvSEFTW.  (47) 


Here,  the  subscript  V  reminds  us  that  we  are  considering  a  variation  at  fixed  volume  (in 
contrast  with  the  original  perturbation  in  our  underlying  quantum  theory).  The  unperturbed 
state  is  taken  to  be  an  equilibrium  vacuum  state.  The  total  entropy  is  therefore  at  an 
extremal  point,  £yStotai  =  0.  A  decrease  in  the  geometric  entropy  (represented  by  the 
boundary  area)  is  thus  compensated  by  an  increase  in  the  entropy  of  the  EFT  state.  Since 
our  original  perturbation  5S  decreases  the  boundary  area  of  our  region,  we  expect  the  field- 


theory  entropy  to  increase.  This  accounts  for  the  minus  sign  in  (46) 


Plugging  (46)  into  (36)  produces  a  relation  between  the  local  scalar  curvature  of  a  region 


around  p  and  the  change  in  the  entropy  of  the  EFT  state  on  the  background: 

2  D(D  +  2  )k\  ££,eft(a) 


anD_1(7<’RC)D+1 


(48) 


Here,  <5S'^T('A^  is  interpreted  as  the  change  in  the  entropy  of  the  EFT  state  in  the  region 
defined  by  Ap,  due  to  shifts  in  both  the  background  geometry  and  the  fields  themselves. 
Considering  nonlocal  rather  than  local  perturbations  results  in  multiplying  the  right-hand 
side  by  i^(0)/2. 


B.  Energy  and  Gravity 


We  can  now  use  the  Entanglement  First  Law  to  relate  the  change  in  entropy  to  an 
energy  density.  The  EFL,  which  relates  changes  in  entropy  under  small  perturbations  to  the 
system’s  modular  Hamiltonian,  holds  true  for  general  quantum  systems.  Given  a  density 
matrix  er,  we  define  its  associated  modular  Hamiltonian  K(a)  through  the  relation 


a  = 


e 


—K 


Tv{e-K) ' 


(49) 


The  Kullback-Leibler  divergence,  or  relative  entropy,  between  two  density  matrices  p  and  a 
is  given  by 


D(p\\a)  =  Tr(plnp)  —  Tr(p  In  cr) 

=  -AS  +  A(K(a)),  (50) 

where 

AS  =  S(p)  -  S{a),  A (K(a))  =  Tr[pK{a)}  -  Tr[aK(a)\.  (51) 

The  relative  entropy  is  nonnegative,  and  is  only  zero  when  the  states  are  identical.  Hence, 
for  infinitesimal  perturbations  cr  =  p  +  dp ,  we  have  D(a\\p)  =  0  to  linear  order,  and  we 
arrive  at  the  EFL  [3j , 

SS  —  S{K).  (52) 
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This  equation  allows  us  to  establish  a  relationship  between  (modular)  energy  and  the  change 
in  entropy,  and  thereby  geometry,  of  our  emergent  space. 

Comparing  to  the  curvature-entropy  relations  for  local  perturbations  (36),  we  see  that 
the  (positive)  induced  local  curvature  Tlp  is  proportional  to  — 5(Kap )■  We  therefore  define 
an  effective  modular  energy  density, 


ep  =  -6(KAp). 


The  curvature  is  then  related  to  the  effective  modular  energy  via 


(53) 


Tip  C  ^Sa  C^pi 


(54) 


where 


2  D(D  +  2) 


«^d-i(7^Rc)£,+1 

is  a  positive  constant.  On  the  other  hand,  an  entangling  operation  in  the  vicinity  of  Ap 
would  give  rise  to  a  negative  “energy”  in  the  region,  adding  negative  spatial  curvature  to 
the  emergent  geometry.  Having  a  large  amount  of  negative  energy  seems  unphysical;  if  our 
model  for  emergent  space  is  to  be  consistent  with  gravity  as  we  know  it,  there  must  be 
conditions  limiting  such  effects,  such  as  instabilities  or  other  dynamical  processes  rendering 
them  unattainable. 


For  nonlocal  perturbations,  we  saw  from  (43)  that  only  positive  curvature  is  generated, 


regardless  of  the  precise  form  of  the  perturbation.  The  curvature-modular-energy  relation 
for  nonlocal  perturbations  is  therefore 


7lp 


(56) 


where  C  =  m(0)|C/2. 

Thus,  once  we  define  an  emergent  geometry  using  mutual  information,  we  see  that  per¬ 
turbing  the  modular  energy  induces  scalar  curvature  in  the  surrounding  space.  This  is 
manifestly  reminiscent  of  the  presence  of  mass-energy  in  the  region  for  the  usual  case  in 
Einstein  gravity.  This  relationship  between  energy  and  curvature  did  not  come  about  by 
“quantizing  gravity” ;  rather,  it  is  a  natural  consequence  of  defining  the  emergent  geometry 
in  terms  of  entanglement. 

However,  the  effective  modular  energy  is  only  an  analogous  expression  for  the  actual 
mass/energy.  To  connect  with  our  familiar  notion  of  energy,  we  need  to  find  the  explicit 
expression  of  the  modular  Hamiltonian  in  terms  of  a  stress  tensor.  Such  expression  will 
generally  be  highly  nonlocal,  and  is  not  known  explicitly  except  for  a  few  cases  [69].  One 
exception,  however,  is  for  a  conformal  held  theory,  where  K  can  be  directly  related  to  the 
stress-energy  tensor  j3j|3].  Consider  a  small  region  centered  at  p  of  size  7^rc,  in  which 

T^FT  is  approximately  constant.  Then  we  have 

x/rr cft\  _  27rHD_ i(7^rc)D+1  x/t,CFT^^\ 

ap  )  —  £)(£)  +  2)  ^00  P  '  57 


Here  5(T^T)  =  Tr[Sp^n>T^J]  is  the  expectation  with  respect  to  the  perturbed  state  of  the 
IR  effective  theory.  Of  course  the  00  component  of  the  stress  tensor  is  simply  the  energy 
density  of  the  theory. 
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Suppose  that  the  RG  flow  JR  for  our  EFT  passes  through  an  IR  fixed  point  at  a  scale  A*. 
At  that  scale,  the  entanglement  structure  of  the  IR  state  can  be  approximated  by  the  ground 


state  of  a  conformal  field  theory,  and  (57)  applies.  Combining  (52)  and  (57)  with  (48)  in 


the  context  of  of  the  low-energy  effective  field  theory,  we  find  that  the  spatial  curvature  is 
related  to  the  EFT  stress  tensor  by 


k  =  E^<:z£ft>, 


a 


(58) 


where  a  is  the  constant  in  the  entanglement /area  relation  (16)  and  K\r  relates  the  full 
entropy  perturbation  to  that  of  the  EFT  as  in  (46).  For  nonlocal  perturbations,  we  saw 


from  (43)  that  only  positive  curvature  is  generated,  regardless  of  the  precise  form  of  the 
perturbation.  The  curvature-modular-energy  relation  for  nonlocal  perturbations  is  therefore 


R  =  ^aWK^cfT) 

( y. 

All  of  our  work  thus  far  has  been  purely  in  the  context  of  space,  rather  than  spacetime. 
We  have  not  posited  any  form  of  Hamiltonian  or  time  evolution.  Let  us  (somewhat  optimisti¬ 
cally)  assume  that  the  present  framework  can  be  adapted  to  a  situation  with  conventional 
time  evolution,  and  furthermore  that  the  dynamics  are  such  that  an  approximate  notion 
of  local  Lorentz  invariance  holds.  (Such  an  assumption  is  highly  nontrivial;  see  e.g.  [70].) 
Thinking  of  our  emergent  space  as  some  spatial  slice  of  a  Lorentzian  spacetime  manifold,  the 
spatial  curvature  can  be  related  to  the  usual  quantities  in  the  Einstein  tensor.  In  particular, 
for  slices  with  vanishing  extrinsic  curvature  we  have 


1ZP  —  2Gqo(p)- 


(60) 


Comparing  to  (58),  we  therefore  find 


Goo  = 


2  ttka* 


6(T& D- 


a 


(61) 


If  we  make  the  identification  27t/va*  / «  — >  8irG,  this  is  nothing  but  the  00-component  of 
the  semiclassical  Einstein  equation.  The  interaction  strength  is  determined  in  part  by  the 
dimensionful  constant  a  that  relates  entropy  to  area,  similar  to  [23].  If  this  reasoning  is 
approximately  true  for  all  time-like  observers  traveling  along  wR  then  one  can  covariantize 
and  arrive  at  the  full  equation 


Gr  =  f1)- 


a 


(62) 


This  is  the  result  for  local  perturbations  using  the  smooth-manifold  approach  to  the  emergent 
geometry.  An  additional  factor  of  F/(0) /2  would  accompany  nonlocal  perturbations  using 
the  lossy  channel  on  Ap,  and  an  analogous  equation  can  be  derived  from  the  Regge  calculus 


approach  using  (B7). 


An  immediate  consequence  of  our  bulk  emergent  gravity  program  is  that  there  is  a  bound 
on  the  change  in  entropy  within  a  region,  reminiscent  of  the  Bekenstcin  and  holographic 
bounds.  We  have  argued  that  positive  energy  corresponds  to  a  decrease  in  the  (full,  UV) 


26 


entropy,  so  we  expect  there  to  be  an  upper  limit  on  the  amount  by  which  the  entropy  can 
decrease.  This  is  of  course  automatic,  as  the  entropy  is  a  nonnegative  number.  Once  the 
region  Ap  is  fixed,  the  maximum  possible  decrease  in  entropy  that  corresponds  to  posi¬ 
tive  “mass-energy”  has  to  be  bounded  by  the  total  entanglement  entropy  of  SAp,  which  is 
proportional  to  the  area  Ap  of  the  region.  More  explicitly, 

I A  SAv\<aAp.  (63) 

This  resembles  the  holographic  entropy  bound.  For  an  entropy  change  that  saturates  the 
bound,  the  vertices  in  regions  Ap  and  Ap  become  disconnected  from  the  graph  point  of 
view.  An  embedding  space  for  Ap  that  reflects  this  change  now  has  a  hole  around  region 
Ap.  Perturbations  that  increase  the  entropy  are  also  bounded,  but  the  bound  scales  with 
the  volume  of  the  region;  we  believe  that  configurations  that  saturate  such  a  bound  do  not 
have  a  simple  geometric  interpretation. 


VI.  DISCUSSION 

We  have  examined  how  space  can  emerge  from  an  abstract  quantum  state  in  Hilbert 
space,  and  how  something  like  Einstein’s  equation  (in  the  form  of  a  relationship  between 
curvature  and  energy)  is  a  natural  consequence  of  this  bulk  emergent  gravity  program. 

We  considered  a  particular  family  of  quantum  states,  those  that  are  “redundancy- 
constrained”  in  a  given  decomposition  of  Hilbert  space.  For  such  states,  a  weighted  graph 
that  captures  the  entanglement  structure  can  be  constructed  from  the  mutual  informa¬ 
tion  between  different  factors,  and  a  manifold  on  which  the  graph  can  be  (approximately) 
isometrically  embedded  is  defined  to  be  its  emergent  geometry.  We  presented  specific 
implementations  of  the  reconstruction  framework  using  the  classical  multidimensional  scal¬ 
ing  algorithm  for  certain  known  area-law  states.  Both  the  dimension  and  the  embedding 
coordinates  for  flat  geometries  can  be  found  through  the  procedure.  At  leading  order,  en¬ 
tanglement  perturbations  backreact  on  the  emergent  geometry,  and  allow  modular  energy 
to  be  associated  with  the  spatial  curvature.  This  relation  is  analogous  to  the  semiclassical 
Einstein  equation. 

A  crucial  feature  of  this  approach  is  that  we  work  directly  with  quantum  states,  rather 
than  by  quantizing  classical  degrees  of  freedom.  No  semiclassical  background  or  asymptotic 
boundary  conditions  are  assumed,  and  the  theory  is  manifestly  finite  (since  regions  of  space 
are  associated  with  finite-dimensional  factors  of  Hilbert  space).  There  is  clearly  a  relation 
with  approaches  that  derive  geometry  from  the  entanglement  structure  of  a  boundary  dual 
theory,  but  the  entanglement  we  examine  is  directly  related  to  degrees  of  freedom  in  the 
emergent  bulk  spacetime.  Because  lengths  and  other  geometric  quantities  are  determined  by 
entanglement,  a  connection  between  perturbations  of  the  quantum  state  and  perturbations 
of  the  geometry  appears  automatically;  in  this  sense,  gravity  appears  to  arise  from  quantum 
mechanics  in  a  natural  way. 

Clearly,  the  framework  is  still  very  incomplete,  and  leaves  much  for  future  investiga¬ 
tion.  An  important  step  in  our  procedure  was  assuming  that  we  were  given  a  preferred 
Hilbert-space  decomposition  "H  =  0  ultimately  we  would  like  to  be  able  to  derive  that 
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decomposition  rather  than  posit  it.  Perhaps  most  importantly,  our  definition  of  distance  in 
terms  of  mutual  information  is  compatible  with  the  behavior  of  field  theories  at  low  ener¬ 
gies,  but  we  would  like  to  verify  that  this  really  is  the  “distance”  we  conventionally  refer 
to  in  quantum  field  theory.  Ultimately  that  will  require  an  investigation  of  the  dynamics 
of  these  states.  An  obvious  next  step  is  to  define  time  evolution,  either  through  the  choice 
of  an  explicit  Hamiltonian  or  by  letting  time  itself  emerge  from  the  quantum  state.  One 
important  challenge  will  be  to  see  whether  approximately  Lorentz- invariant  dynamics  can 
be  recovered  at  low  energies,  and  whether  or  not  the  finite  nature  of  Hilbert  space  predicts 
testable  deviations  from  exact  Lorentz  symmetry.  We  might  imagine  that,  given  a  state  \i(j) 
whose  geometry  is  constructed  using  entanglement,  one  can  generate  all  time-slices  using  a 
known  local  Hamiltonian  such  that  \i[})  is  a  low  energy  state.  Alternatively,  by  working  with 
mixed  states  one  could  adopt  the  thermal  time  hypothesis  in]  and  generate  state-dependent 
time  flow  purely  from  the  modular  Hamiltonian,  which  is  in  principle  attainable  from  just 
the  density  operator. 

To  analyze  the  emergent  geometries  of  states  beyond  redundancy-constraint,  deeper  un¬ 
derstandings  of  the  entropy  data  for  subregions  of  different  sizes  will  be  important.  One  such 
case  is  manifest  in  the  context  of  AdS/CFT  correspondence,  where  entanglement  entropy 
of  different-sized  balls  in  the  CFT  are  needed  to  obtain  bulk  geometric  information  through 
a  radon  transform.  One  such  approach  may  be  to  introduce  additional  structures  on  the 
graph  and  extend  it  to  a  tensor  network.  The  program  of  geometry  from  tensor  networks 
has  mostly  been  based  on  states  with  a  high  degree  of  symmetry,  such  that  notions  of  length 
and  curvature  can  be  assigned  through  simple  geodesic  matching  and  tessellation  of  space. 
The  results  obtained  here  suggest  that  for  tensor  networks  with  small  perturbations,  one 
can  also  modify  the  geometric  assignment  accordingly,  matching  the  change  in  correlation 
or  entanglement  to  perturbation  in  geodesic  lengths.  A  notion  of  (coarse)  local  curvature 
can  also  be  defined  on  triangulated  spaces  using  entanglement  and  Regge  calculus,  which 
seems  more  natural  for  programs  that  relates  network  geometries  to  those  of  spacetime. 

The  emergence  of  time  evolution  will  also  be  useful  for  the  study  of  more  complex  be¬ 
haviors  related  to  entanglement  perturbations.  For  instance,  one  can  examine  the  inter¬ 
actions  among  multiple  perturbations  created  in  some  local  region.  If  the  model  is  truly 
gravitational,  the  time  evolution  experiment  should  be  consistent  with  our  knowledge  of 
gravitational  dynamics.  It  will  also  be  interesting  to  study  the  redundancy-constrained  de¬ 
formations  of  states  beyond  perturbative  limit.  Intuitively,  we  expect  the  emergence  of  a 
classical  wormhole  geometry  by  nonlocally  entangling  large  number  of  degrees  of  freedom  in 
a  coherent  manner.  One  can  also  examine  purely  quantum  phenomena  outside  the  context 
of  classical  Einstein  gravity,  including  black-hole  entropy  and  evaporation,  using  mutual 
information  rather  than  classical  spacetime  geometry. 
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Appendix  A:  Redundancy-Constraint  and  Coarse-graining 


In  this  appendix  we  consider  how  to  construct  a  coarse-grained  decomposition  of  Hilbert 


space  that  is  redundancy-constrained  (RC),  as  defined  in  Section  III  A  when  an  initial  fine¬ 
grained  one  is  not. 

Given  some  state  |Ao)  £  77  and  some  fixed  Hilbert  space  decomposition  77  =  77,;  for 

M  sufficiently  large,  we  first  create  a  network  represented  by  a  graph  G0  =  (V0,E0).  The 
graph  has  N  vertices  labelled  by  i,  and  each  edge  {7,  j}  is  weighted  by  I{i  :j),  where  I(i  :j) 
is  the  mutual  information  of  partitions  i  and  j.  If  the  resulting  graph  is  RC  to  the  desired 
degree  of  accuracy,  no  further  coarse-graining  is  needed. 

If  not,  consider  the  set  of  all  partitioning  schemes  C  for  a  coarse-grained  decomposition 


of  the  Hilbert  space  such  that  77  =  (££)*  77p.  For  each  scheme  S  6  C,  we  require  that  N  <  M 
for  some  sufficiently  large  N  so  that  non-trivial  entanglement  structure  is  still  allowed.  Each 
partition  S  =  {{A,  A,  •  •  •  },  {A,  A+i, corresponds  to  constructing  a  more  coarse¬ 
grained  decomposition  of  the  Hilbert  space  by  taking  the  union  of  original  subfactors;  that 
is,  for  each  sp  G  S,  77p  =  <S)ieSp  'Hi- 

A  partition  is  RC-valid  if  the  mutual-information-weighted-network  G  =  (V,  E)  based 
on  the  coarse-grained  decomposition  is  redundancy-constrained.  While  there  is  no  obvious 
way  to  choose  the  best  coarse-graining  scheme  at  this  point,  it  is  natural  to  consider  the 
most  uniform  partitioning,  so  that  all  Hilbert-space  subfactors  have  approximately  equal 
dimensions. 

If  no  such  coarse-graining  can  be  found  for  N  reasonably  large,  then  the  procedure  fails 
and  we  are  forced  to  conclude  that  the  given  state  cannot  be  cast  in  RC  form  in  a  simple 
way.  We  do  not  claim,  however,  that  it  doesn’t  admit  a  simple  geometric  description,  as  this 
is  clearly  false  from  our  knowledge  of  AdS/CFT.  Reconstructing  the  geometry  from  such 
states  are  interesting  problems. 

A  more  specific  method  for  coarse-graining  can  also  be  constructed  using  network  renor¬ 
malization.  While  the  algorithm  is  less  computationally  intensive,  it  fails  for  certain  states 
if  the  original  decomposition  yields  little  useful  information.  For  instance,  it  fails  for  the 
ground  state  of  Toric  code  on  a  square  lattice  where  the  given  Hilbert  space  decomposition  is 
the  usual  tensor  product  of  spin-1/2  degree  of  freedom  on  each  link.  It  is,  however,  useful  for 
certain  finitely  correlated  states  and/or  those  of  typical  condensed  matter  system  at  scales 
larger  than  the  correlation  length.  It  also  works  for  the  toric  code  if  the  decomposition  is 
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more  coarse-grained. 

For  such  states  and  decompositions,  we  follow  network  renormalization  procedure  by 
again  constructing  the  mutual-information-weighted-network  Go  =  (Vo,  Do)-  Assume  Go  is 
connected;  in  the  case  when  Go  has  multiple  large  disconnected  components,  one  can  simply 
perform  the  procedure  separately  for  each  connected  component. 

Then  proceed  to  define  a  metric  do(i,j)  on  Go  in  the  same  manner  as 
Then,  for  any  vertex  v  on  the  graph,  we  consider  an  e-ball  Be(v )  such  that 


14),  with  £rc  -3  £q. 


Be(v)  =  {v'  G  V\ d(v,v')  <  e}, 


(Al) 


where  d(v,v')  is  a  metric  defined  on  the  set  of  vertices  V.  Seed  the  entire  graph  with 
points  like  v  until  the  whole  graph  is  covered  by  e-balls.  Choose  a  minimum  cover  and 
compute  the  entanglement  entropy  for  each  the  union  of  subregions  in  each  ball  to  generate 
a  coarse-grained  graph  G i,  where  each  vertex  now  is  labeled  by  the  union  of  the  original 
subregions  and  the  edge  weights  are  given  by  the  mutual  informations  of  the  coarse-grained 
subregions.  Repeating  the  coarse-graining  procedure  until  all  higher  subregions  B\  can  be 
well  approximated  by  the  cut  function  (J9]) ,  we  then  label  the  coarse-grained  subregion  at 
this  scale  to  be  Ap  for  p  G  S,  and  label  the  corresponding  coarse-grained  network  with  Ap 
also.  The  resulting  state  will  be,  to  a  good  approximation,  redundancy-constrained. 


Appendix  B:  Entanglement  Perturbations  and  Coarse  Curvature 


In  this  appendix  we  consider  an  alternative  approach  to  calculating  the  curvature  induced 
by  an  entanglement  perturbation,  working  directly  with  the  discrete  graph  rather  than 
finding  an  embedding  Riemannian  manifold.  Here  we  use  the  techniques  in  Regge  calculus 
in  which  the  sense  of  spatial  curvature  is  determined  by  deficit  angles. 


For  a  space  of  fixed  integral  Hausdorff  dimension  D  obtained  from  (15),  consider  a  vertex 
p  and  construct  a  local  triangulation,  if  one  exists.  We  say  the  space  is  r-locally  triangulable 
at  p  if  one  can  construct  a  abstract  simplicial  complex  /C,  where  the  simplices  are  sets  of 
vertices  in  the  metric  subspace  near  p,  by  imposing  a  distance  cutoff  r,  and  if  there  exists 
an  isometry  (with  respect  to  the  metric  distance)  that  maps  /C  to  a  geometric  simplicial 
complex  K  where  inter-vertex  distances  of  the  metric  subspace  are  preserved.  If  K  is  also  a 
simplicial  manifold  we  can  proceed  to  define  Regge  curvature. 

Select  a  codimension-2  simplex  X  3  p  as  a  hinge,  its  volume  given  by  (18).  As  the 


simplices  in  K  are  equipped  with  the  usual  Euclidean  inner  product,  angles  can  be  defined 
and  deficit  angle  at  the  hinge  8 (A)  =  27T  —  0(X)  can  computed  using  the  inner  product 
structure  in  Euclidean  space.  Here  we  define 


»P0  =  X>P0. 


(Bl) 


where  4>i(X)  is  the  angle  between  the  unique  two  faces  of  a  simplex  containing  the  hinge  X. 
In  the  case  of  a  D-dimensional  area-law  system  where  I(p:q )  ^  0,  construct  a  simplex  by 
considering  the  n  shortest  distances.  The  curvature  is  then  related  to  the  deficit  angles  5* 
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in  the  region  by 


Rt  =  E  S*L‘ • 


(B2) 


i 

where  are  the  volumes  of  the  codimcnsion-2  hinges  at  which  the  curvature  is  concentrated. 
For  D  =  2,  the  hinge  is  a  point,  and  the  curvature  is  given  by  the  deficit  angle,  where  we  set 
Li  =  1.  In  the  continuum  limit,  (|B2[)  becomes  f  dDx^/g1Z,  where  7 7  is  the  scalar  curvature. 
In  the  case  of  emergent  Euclidean  (flat)  space,  we  require  Rt  =  0. 

Let  us  consider  the  effect  of  an  entanglement  perturbation  on  the  geometry  of  a  distance- 
function  graph  G,  using  this  technique.  Again,  the  forms  for  both  local  and  nonlocal  per¬ 
turbations  are  similar  up  to  an  overall  factor. 

Since  the  original  deficit  angle  Sp  =  0  in  flat  space,  for  p  lying  on  the  hinge,  we  may 
consider  the  angular  deficit  produced  by  varying  each  elongated  edge  connected  to  p.  For 
each  such  simplex  S  connected  to  the  hinge,  the  deficit  at  p  induced  by  varying  the  length 
lj  of  each  edge  in  the  simplex  assumes,  at  leading  order,  the  form 

1  =  — t  +  om/ifn  <B3) 

where  lj  denotes  the  jth  edge  length  of  the  /^-simplex  S.  t \S\h ,  h,  ■  ■  ■ )  is  a  function  that 
has  dimension  of  length  and  depends  on  the  edge  that  is  varied,  as  well  as  all  the  edge 
lengths  that  connect  the  vertices  of  S. 

The  overall  deficit  in  a  triangulation  where  all  edges  have  roughly  the  same  (unperturbed) 
length  £  is 

= £  4  +  omnf?)  (B4) 

j,s 

=  Ur(D)j  +  0((5(/eYj,  (B5) 


where  A fp(D)  depends  on  the  simplices  in  the  triangulation.  For  example,  if  equilateral 
triangles  with  sides  £  are  used  to  triangulate  the  2-dimensional  flat  space  around  p,  then 
A fp(D  =  2)  =  12/ \/3  in  the  case  where  all  edges  emanating  from  p  vary  by  the  same  amount 
8£  under  the  entanglement  perturbation.  Note  that  in  dimension  greater  than  2,  there  is  no 
uniform  tiling  such  that  all  edges  are  equal,  hence  the  approximation  in  some  average  sense. 

To  estimate  the  coarse  curvature  at  p,  take  Lp  ~  £D~2  as  the  volume  of  the  codimension- 
2  hinge.  S£  is  identified  with  5d(p,  q)  for  the  change  in  distance  between  adjacent  vertices, 
and  £  ~  d(p,  q)  for  all  q  immediately  adjacent  to  p  in  the  triangulation.  The  total  coarse 
curvature  is  given  by 

Rc  =  A PLP  ~  A rp(D)£D~35£.  (B6) 


On  the  other  hand,  we  know  from  [72]  that  in  the  continuum  limit,  if  the  metric  g  is 
approximately  constant  in  the  small  region,  Rc  — »  £D1Z,  where  1Z  is  the  average  coarse 
scalar  curvature  of  the  space  contained  in  the  small  region  with  approximate  volume  iD . 
Applying  (39)  and  the  EFL,  the  average  coarse  curvature  in  the  region  is 


77  =  ZS{H). 


(B7) 
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Here,  for  a  perturbation  induced  by  a  unitary  UApj  ,  the  constant  Z  =  £yiq£(p,  Q-  D)  depends 
on  the  triangulation,  strength  of  entanglement,  and  choice  of  coarse-graining.  The  factor 
£(p,  q,  D )  oc  \Afp(D)&'(i(p :  q))\/I0(p :  q),  which  parametrizes  all  the  order-one  constants  that 
enter  into  the  process  of  averaging,  can  be  explicitly  computed  once  the  triangulation  and 
the  inverse  function  $  are  known.  We  have  taken  l  =  yfac  to  denote  the  average  radius 
of  the  region,  as  before.  Note  that  this  is  consistent  with  (40)  up  to  the  dimension-  and 
triangulation-dependent  factors. 


For  a  nonlocal  perturbation  through  a  channel  A,  because  the  unitary  only  acts  on  Ap 
and  the  ancilla,  the  mutual  information  Si (s  :  q)  =  0,  and  hence  the  distance  functions  d(s,  q) 
are  invariant  for  all  s,q  7^  p.  Only  legs  emanating  from  the  vertex  p  in  the  triangulation  are 
varied.  The  unitary  perturbing  map  has  a  wider  range  of  possible  consequences,  as  its  form 
is  unspecified.  Although  the  values  of  Z  may  be  different  depending  on  the  specific  map 
used,  the  formalism  remain  the  same.  As  a  result,  for  nonlocal  perturbations  with  channel 
A,  we  have  an  equation  of  the  same  form,  with  Z  — >  ZF'A( 0)/2.  In  particular,  if  we  restrict 
U apap  only  remove  entanglement  symmetrically  near  the  boundary  of  the  region  Ap  as 
before,  the  values  of  Z  will  be  the  same. 
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